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Abstract 


We consider an XY quantum spin chain that consists of a left, center and right part initially 
at thermal equilibrium at temperatures Ti, T c , and T r , respectively. The left and right systems 
are infinitely extended thermal reservoirs and the central system is a small quantum system linking 
these two reservoirs. If there is a temperature differential, then heat and entropy will flow from 
one part of the chain to the other. We consider the Evans-Searles and Gallavotti-Cohen functionals 
which describe the fluctuations of this flux with respect to the initial state of the system and the 
non-equilibrium steady state reached by the system in the large time limit. We also define the 
full counting statistics for the XY chain and consider the associated entropic functional, as well a 
natural class of functionals that interpolate between the full counting statistics functional and the 
direct quantization of the variational characterization of the Evans-Searles functional which appears 
in classical non-equilibrium statistical mechanics. The Jordan-Wigner transformation associates a 
free Fermi gas and Jacobi matrix to our XY chain. Using this representation we are able to compute 
the entropic functionals in the large time limit in terms of the scattering data of the underlying Jacobi 
matrix. We show that the Gallavotti-Cohen and Evans-Searles functionals are identical in this limit. 
Furthermore, we show that all of these entropic functionals are equal in the large time limit if and 
only if the underlying Jacobi matrix is reflectionless. 
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Preface 

This document is a thesis submitted in partial requirements of the degree of M. Sc. in Mathematics, 
sought by the author, Benjamin Landon. The author was supervised by Dr. Vojkan Jaksic and Dr. 
Robert Seiringer. The thesis has been written by the author. Dr. Jaksic and Dr. Claude-Alain Pillet 
provided guidance where necessary as to the content of the thesis and also suggested methods by which 
to complete the proofs. 

1 Introduction 

Entropic fluctuations in statistical mechanics concerns the study of the relative probability that entropy 
will increase or decrease over a time interval in a physical system. We provide a brief historical review 
of the subject; the interested reader may consult [RM] for a more exhaustive treatment and additional 
references. The study of fluctuations in statistical mechanics began in 1905 with the seminal paper of 
Einstein [El] on Brownian motion in which the first fluctuation-dissipation relation was given. In his 
1910 paper [E2], Einstein gave a link between the entropy formula of Boltzmann and the probability of 
fluctuations out of an equilibrium state. Subsequent developments were made by Ornstein and Nyquist; 
Ornstein obtained a fluctuation-dissipation relation for a random force acting on a Brownian particle [Or] 
and Nyquist computed spectral densities and correlation functions of thermal noise in linear electrical 
circuits in terms of their impedance [Ny]. The classic result of Onsager is the converse [01, 02]; he 
obtained a formula for the transport coefficients or impedances in terms of thermal fluctuations. Progress 
in the theory of transport coefficients and fluctuation-dissipation relations was continued in the works of 
Green [Gl, G2] and Kubo [Kub] (e.g., the derivation of the well-known Green-Kubo formulas). 

These early results concerned primarily the close-to-equilibrium regime; e.g., the transport coefficients 
refer to the first order response of a physical system in equilibrium to thermodynamic or mechanical forces 
pushing it out of equilibrium. While the equilibrium theory is considered satisfactory and complete, the 
same is not true of the non-equilibrium theory [RM] . 

Modern fluctuation theorems were first suggested numerically in the paper [ECM] . There, the authors 
were interested in violations of the second law of thermodynamics; in the course of studying a deterministic 
particle system they found the relationship (roughly speaking) that 

W) = e-VPti-fi (1-1) 
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where Pt{4>) is the probability of measuring a decrease in the entropy of the system of magnitude <j> over 
a time interval of length t. That is, decreases of entropy can occur, but the probability of an entropy 
decrease occuring is exponentially small in the magnitude of the decrease compared to measuring an 
increase of the same magnitude. This exponential suppression has an equivalent formulation in terms 
of the cumulant generating function E t (a ) for the random variable describing the entropy change over 
the time interval [0, £]. To be more precise, (1.1) holds iff the cumulant generating function obeys 
E t (a) = E t { 1 — a); such a relationship has come to be known in the literature as a fluctuation relation. 
Such fluctuation relations were established theoretically for the first time by Evans and Searles [ IS 
and Gallavotti and Cohen [GC1, GC2]. One of the novel features of such modern fluctuation theorems 
is that they hold for systems far from equilibrium, and in the linear regime near equilibrium reduce 
to Green-Kubo formulas and Onsager relations, an observation made by Gallavotti in 1996 [Ga, JPR], 
They are therefore a generalization or a far-from-equilibrium version of the Green-Kubo formulas and 
Onsager relations which have played an important role in the development of non-equilibrium statistical 
mechanics. 

While the majority of work in entropic fluctuations concerns classical mechanics, comparatively little 
is known in the quantum case and establishing quantum analogues of the existing results is an emerging 
and active area of research. The modern fluctuations theorems represent some of the few exact and general 
results in non-equilibrium statistical mechanics and it is therefore desireable to obtain a generalization 
to the quantum regime. As it stands, both the classical and quantum theory admit an axiomatization 
[JPR, JOP]; starting with a general classical or quantum dynamical system, the basic objects of the 
theory such as the entropy production observable and the finite time entropic functionals can be defined 
at a great level of generality. The axioms concern the existence and regularity of the entropic functionals 
in the large time limit. In the quantum case, it is generally extremely difficult to verify these axioms in 
physically interesting models [JOPP]. 

The mathematical tools required to deal with the quantum case are in general very demanding. Virtu¬ 
ally all aspects of Tonrita-Takesaki modular theory play a role in the description of entropic fluctuations 
in quantum non-equilibrium statistical mechanics. For example, the Araki-Masuda non-commutative L p 
spaces take the place of the familiar L p spaces of measure theory appearing in the classical case, and the 
Connes cocyclc and relative modular operators take the place of Radon-Nikodyn derivatives [JOPP]. As 
the quantum case is somewhat overshadowed by these technical requirements, it makes sense to consider 
a specific, exactly solvable model, such as the XY spin chain, and compute the quantities of interest in 
this specific case. Due to the exact solvability, the complete description of the entropic fluctuations of 
the XY quantum spin chain will require no more technical tools than the basics of trace-class scattering 
theory, a subject accessible to beginning graduate students or even advanced undergraduates. The XY 
spin chain has seen a great deal of study in quantum statistical mechanics - we refer the reader to [JLP] 
for a sampling of the existing literature on the subject. For example, the first proofs of the existence of 
a non-equilibrium steady state were given in the context of XY chains [AH, AP]. 

As previously mentioned, the goal of this thesis is to test the emerging theory of entropic fluctuations 
in the context of XY quantum spin chains. In particular, we compute the large time limit entropic 
functionals of the XY chain. These functionals describe the fluctuations of entropy production in the XY 
spin chain. The results of this paper were already announced in [JLP] , and some proof sketches appeared 
there. The model considered here is more general, and we complete the calculations that were omitted 
in [JLP]. Our principle results concern the existence and the properties of the entropic functionals and 
are the content of Theorem 4.10; we summarize the results in the sequel. 

1.1 Outline of paper 

This paper is organized as follows. Section 2 deals with the finite volume XY chain and studies some of 
its properties. We initially define an XY chain confined to a finite interval in Z. We then define the finite 
volume open XY chain which consists of three pieces, each of which is by itself an XY chain confined 
to a finite interval. We refer to these pieces as the left, center, and right parts of the chain. The parts 
are coupled at their end points; initially the left/center/right part of the chain is at thermal equilibrium 
at inverse temperature /3 i/ c / r . If there is a temperature differential, then energy and entropy will flow 
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from one part of the chain to the other. We study the properties of this flux and associated entropy 
production, and introduce the finite time and finite volume entropic functionals. 

The entropic functionals introduced are the following. The Evans-Searles functional is the direct 
quantization of the corresponding functional in classical non-equilibrium statistical mechanics. In the 
classical case, the functional describes fluctuations of the entropy production of the system with respect 
to its initial state. However, the direct quantization results in a functional that fails to have the desired 
a <-> 1 — a symmetry characterizing the aforementioned fluctuation relation of non-equilibrium statistical 
mechanics (see Section 2.3.2). 

At this point we introduce an entropic functional associated to the full counting statistics of the XY 
chain, which is a probability measure associated to a repeated quantum measurement protocol of the 
entropy of the XY chain (see Section 2.3.3 for more details). We introduce a functional which is the 
direct quantization of the variational characterization of the Evans-Searles functional in classical non¬ 
equilibrium statistical mechanics (see Proposition 2.3(vii)), and a class of functionals which interpolate 
between the variational functional and the full counting statistics functional. We mention that the full 
counting statistics, variational and interpolating functionals all have the a <-» 1 — a symmetry. 

In the latter half of Section 2 we introduce the Jordan-Wigner transformation. The Jordan-Wigner 
transformation maps the confined XY chain to a free Fermi gas, and associates to the Hamiltonian of 
the XY chain a Jacobi matrix. In the remainder of the paper, we work exclusively in the Fermi gas 
representation of the XY chain. The Jordan-Wigner representation allows us to derive simple formulas 
for the entropic functionals which prove useful in the remainder of the paper. 

In Section 3 we take the thermodynamic limit in which the volume of the chain tends to infinity. 
In our set-up, the center part of the chain is fixed, while the right and left parts of the chain become 
semi-infinite chains (that is, they have one fixed end, to which the center part is coupled, and one end 
at infinity). We call the resulting object the extended XY chain. We show that the initial state and 
dynamics of the finite volume chains converge to a state and dynamics of the extended XY chain. We 
also show that the finite volume entropic functionals have infinite volume limits and compute them in 
closed form. 

Section 4 concerns the large time limit, and contains our main results, Theorem 4.10. In the previous 
section we fixed the time t and took the thermodynamic limit, resulting in, for each fixed f, an extended 
XY chain. In this section, we will take the limit t —» oo of the extended XY chain. We begin by introducing 
the necessary facts from spectral and scattering theory that we will require to prove our resuts. In the case 
of the XY chain, the large time limit is easy to control via trace-class scattering theory. We prove that in 
the large time limit, the state of the XY chain converges to a non-equilibrium steady state (NESS) and 
we compute the steady state heat fluxes in terms of the scattering data of the underlying Jacobi matrix 
of the XY chain. Except in trivial cases, there is a non-zero steady state entropy production as long as 
the left and right parts of the chain are initially at different temperatures. At this point we introduce 
the Gallavotti-Cohen entropic functional which measures fluctuations with respect to the NESS. 

Subsequently, we show that the entropic functionals have large time limits and compute them in closed 
form in terms of the scattering data of the Jacobi matrix associated to the XY chain. We prove that, 
except in trivial cases, all the entropic functionals are identical in the large time limit iff the Jacobi matrix 
associated to the XY chain is reflectionless (the definition of which is given in Section 4.3). If the Jacobi 
matrix is not reflectionless, then the full counting statistics, variational and interpolating functionals are 
all different, and are different from the Gallavotti-Cohen and Evans-Searles functionals. Regardless of 
whether or not the underlying Jacobi matrix is reflectionless, the Evans-Searles and Gallavotti-Cohen 
functionals are equal in the large time limit. In particular, the Evans-Searles functional recovers the 
a <-> 1 — a symmetry if and only if the Jacobi matrix is reflectionless. 


2 Finite volume XY chain 

2.1 The XY chain confined to an interval 

A finite dimensional quantum system is a pair ( TL , O-h) where T-L is a finite dimensional complex Hilbert 
space and On is the algebra of matrices acting on T-L. On is called the algebra of observables. Throughout 
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this paper the inner product on a Hilbert space T-L will be denoted (•. ■')r and will be taken to be linear 
w.r.t. the second variable. When the meaning is clear we will drop the subscript and write instead (•, •). 
The self-adjoint elements of Or are those A s.t. A = A* and the positive elements of Or are those A 
s.t (ip, Aip ) > 0 for every ip e T~i. The spectrum of A e Or is denoted sp(A) and is the set of complex 
numbers z s.t. A — z is not invertible. 

If A is an eigenvalue of A, we denote the projection onto the corresponding eigenspace by 1 a (A). We 
set | A| = VA*A and for p e]0, oof set 

\\A\\ p =(tr\A\n 1/p . ( 2 . 1 ) 

Furthermore, we define HA^ to be the largest eigenvalue of |A|. It follows that lim^oo M P = Moo- 
We now consider an XY spin chain confined to a finite interval in Z. To each xeZwe associate the 
Hilbert space 1~L X = C 2 . The algebra of matrices acting on T-L x is denoted O x and is equal to M 2 (C), that 
is, 2 x 2 matrices with complex entries. A basis for O x given by the Pauli matrices 


cxW = 


0 1 
1 0 


together with the identity matrix l x . 



The Pauli matrices satisfy the relations 
= 8 jk l x + ie 3kl a^\ 


( 2 . 2 ) 


(2.3) 


with Sjk the Kronecker delta and e 3kl the Levi-Civita symbol. 

Let A = [N, M ] be a finite interval in Z. The XY chain confined to A is the quantum system described 
by the Hilbert space and algebra of observables 

n A = (g)H x , 0 A = (g)0 x . (2.4) 

xeA xeA 

For simplicty of notation, we identify the operator A x e O x at a single site with the operator (<8> y6 A\{x}li/)® 
A x which is an element of O a , the algebra of observables of the confined XY chain. Similarly, if A' c A, 
we identify an element A e 0 A ' with the element A® 1 a\a' e O a . The Hamiltonian of the XY chain 
confined to A is given by 


H a 


xe[N,M[ 


+ oY’o 


(2U 2 ) 


x+1 



(2.5) 


Here, {J x } xe z and {X x } xe z are assumed to be bounded sequences of real numbers. J x is the nearest 
neighbour coupling at each site, and Ax is the strength of a magnetic field in the direction (3) at the site 

x. 


2.2 Finite volume open XY chain 

In the previous section we described an XY chain confined to an interval in Z. The open XY chain is an 
XY chain confined to an interval where the initial state of the system consists of three distinct pieces, 
each at thermal equilibrium at a different temperature. This will be made precise in this section. 

Let A = [— M, M ] be a finite interval in Z, and consider an XY chain confined to A as defined in Section 
2. Let N e Z + (here, Z + = N u {0}), with N < M. We denote = [— M, —N— 1], A c = [— N, X] and 
A r = [X + 1, M]. A l / C / R is the left/center/right part of the chain. The Hamiltonian of the XY chain 
confined to A can be written as 


H A = Hl + He + Hr + V L + Vr (2.6) 

where H l /c/r is the Hamiltonian of an XY chain confined to A l/c/r, as defined by (2.5), and 

Fi = \ j-N -+ CT-N-^-n), Vr = ^Jn^n^N+I + ^N^N+l), ( 2 - 7 ) 
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is the coupling energy of the left and right chains to the center. We define H 0 = Hr + He + Hu and 
V = Vl + Vr. Furthermore, until the end of the section, A = [— M, M ] is fixed and we will omit the 
subscript A when the context is clear (i.e., we write H = H\. H = "Ha, etc.) 

A density matrix p on the algebra of observables of a quantum system is a positive matrix with 
tr(p) = 1. The state associated to the density matrix p is the linear functional on the algebra of 
observables given by p{A) = tr(pA), for A an observable. Since every positive linear functional v acting 
on the algebra of observables with i/(l) = 1 is given by v(A) = tr( ia 4), with v a density matrix, we will 
abuse the notation slightly and identity the state with the density matrix. 

If A is an algebra of observables of a quantum system, then a ^-automorphism of A is a linear bijection 
d : A A so that i?( AB) = , d(A)‘d(B) and d(A*) = 'd(A)*. It follows that sp(d(A)) = sp(A) and that 
||'d(A)|| = ||A|j. The set of ^-automorphisms of A is denoted Aut(_4). 

A dynamics on the algebra of observables A of a quantum system is a continuous one-parameter 
subgroup of ^-automorphisms of A. That is, it is a map IR 3 t —» r* e Aut(_4) so that r t+s = r* o r s and 
lim t ^o It 4 (A) — ^4|| = 0 for every A in A. If B is a self-adjoint element of A, then r*(A) := e ltB Ae~ ltB 
is a dynamics on A. 

The Hamiltonian of the XY chain induces the dynamics on the algebra of observables, defined for 
A e O as 


t\A) = e itH Ae~ itH . (2.8) 

In the Heisenberg picture observables evolve forward in time and in the Schrodinger picture states evolve 
backwards in time and so we define, 


A t = r‘(A), p t = T \p). 

With this convention, pt{A) = p{A t ). 

The initial state of the open XY chain is 

g—/3 lHl-\ — PcHc—PrHr 
W tr{e ~^ LHL ^— PcHc—PrHrY 

With this initial state, the left, center and right systems are initially at thermal equilbrium at inverse 
temperatures /3l,Pc and /3fl, respectively. In fact, in the absence of the coupling between the chains, 
the state w is a steady state. That is, if the dynamics in (2.8) were induced by Hq instead of H, then oj 
would be invariant under the dynamics as [Hq, uj] = 0 (here, [A, B] = AB — BA denotes the commutator 
of observables A and B ). Note also that w is a faithful state. That is, kerw = {0}. We shall call the 
quantum dynamical system defined on the Hilbert space H the open XY chain confined to A. 

From now until the end of the paper we take, 


(2.9) 


( 2 . 10 ) 


Pc = 0. 


( 2 . 11 ) 


This convention is for notational simplicity and does not affect our main results. 

We also note that the XY chain is time reversal invariant (TRI). That is, there exists an anti-linear 
involutive ^-automorphism 0 of O s.t. r* o 0 = 0 o r _t and 0(w) = 0. 0 is called a time reversal and 
will be described in Section 2.4. 

For future reference, we define the relative entropy of a state p with respect to a state v to be 


S(p\v) 


{ tr(p(log v — log p)) if ker v c ker p, 
—oo otherwise. 


( 2 . 12 ) 


With the convention that 0 x oo = 0, the relative entropy is well-defined and as a consequence of Klein’s 
inequality (see Theorem 2.1 in [JOPP]), S(p\v) < tr {v — p) = 0 with equality iff p = v. 
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2.3 Observables of the open XY chain 

2.3.1 Entropy production and heat fluxes 

In this section we define some of the basic observables of the open XY chain confined to an interval. 
Much of the discussion here coincides with the abstract framework put forth in [JOPP] in the case of 
general finite dimensional quantum systems, and our definitions are the same as those in [JLP]. As the 
Hamiltonian Hl/r is associated with the energy in the left/right part of the chain, the heat fluxes out of 
the left and right parts of the chain at t, = 0 are 


®l/r ~ 




t =0 


—i[H, H l / r ] - i [H l /r,V l / r \. 


(2.13) 


It is easy to compute these fluxes using the relations (2.3): 


*7 T / (1) r (2) (l) (2) (2)r 

= -^JnJn + 1 [<?N l a N+l’ a N+li a N+2 + a N P 7 


A 2 ) U 1 ) 


N L U JV+1’ u N + li u N+2 


T (3) <t (2) 

J N + 1> a N+l 


1) 


+ ^ JjvAjv+1 (ajvVw>i’°vli] +< J N ) [° 

-Ir T rr (3) fn-WH 2 ) - J 2 ) J 1 ) \ A- - T \ /A 2 ) A 1 ) _ n -( 1 ) n -( 2 ) 

— r> 'JN’JN+l <J N + l \y a N a N +2 a N a N+2J + ^Jn^N + I ^JV + l a N a N+lJ ■ (2.I4J 


A similar computation shows that 


1 ( 3 ) 

= y u — ]\f—2 U — N 


r (2) 


( 1 ) 

o'! m - cr 


(i) 

— N—2 U — 


fT ~^) + (a 


(2) (1) 

N°—N — l 


(!) _( 2 ) 


Cr _ 


The entropy production observable (as defined by thermodynamics) is 

cr = -Pl®l ~ Pr$r = -i[H,\oguj] 


-N u -N-lj 

(2.15) 


(2.16) 


Note that the heat fluxes and the entropy production observable change sign under the time reversal 0: 

e(*L/fl) = -$L/R, 0 (a) =-a (2.17) 

This is a consequence of the time reversal invariance of u. The mean entropy production rate over the 
time interval [ 0 ,f] is 

1 r 4 

E 4 = - cr s ds. (2.18) 

t Jo 

If t < 0, we define X* by the above equation as well. We summarize a few of the basic properties of these 
observables in the proposition below. These results as well as further details can be found in [JOPP] (see 
also [JLP]). 

Proposition 2.1. The folowing hold: 

(i) logw t = logo; + tT~ t ( E 4 ), 

(ii) S(w t |u;) = -fw(E 4 ), 

(hi) t 4 (X -4 ) = X 4 , 

(iv) E 4 = -r 4 (0(E 4 )), 

(v) The spectrum of X 4 is symmetric with respect to 0. In fact, diml^cd) = diml_ 0 (cr 4 ) for all 
</GS P (X t ). 

Proof, (i) We differentiate and integrate to obtain, 

log uj t - logo; = J ^-logw s ds = J r~ s (-i[H, logo;]) = tr” 4 (E 4 ). (2.19) 
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(ii) Using (i), we have 


5(wt|o;) = wt(logw — logwj) = — ta;(r t (r 4 (£ 4 ))) = —tu (£ 4 ). (2.20) 


(iii) Using (i), we have t 4 (£ 4 ) = (logw — logw_ t )/f. On the other hand, (i) implies that E 4 = T 4 (logw t — 
logw)/f = (log D - log UJ—t)/t. 

(iv) By (2.17), 0(£ 4 ) = —E -4 . By (iii), £ 4 = r 4 (E" 4 ) = r 4 (0(0(£- 4 ))) = -r 4 (0(E 4 )). 

(v) By Exercise 3.1 in [JOPP] (see also [ ]), 0(A) = UqAUq 1 for an anti-unitary Uq. Therefore (iv) 

implies that for every eigenvector v with eigenvalue (j> of E 4 , e ltH Ue>v is an eigenvector of £ 4 with eigenvalue 
—</>, and the claim follows. □ 

Proposition 2.1 has a few consequences. Firstly, part (i) allows one to interpret the entropy production 
observable as the quantum phase space contraction rate. Part (ii) together with the fact that the relative 
entropy of two states is always nonpositive implies that the average entropy production over the interval 
[0,<] is nonnegative. By this fact, 


0 < fw(£ 4 ) 


-PL 


f w($L jS )ds - p R f w($ fljS )ds = -f 3 L A L ,t ~ PrArj 
Jo Jo 


( 2 . 21 ) 


where A L / Rt = Jq w($ i / i?)S )ds is the average heat flow out of the left/right part of the chain during the 
time interval [0, t]. This inequality is the finite time second law of thermodynamics; on average, heat 
flows from the hotter to the colder part of the chain. 

Additionally, parts (ii) and (v) of Proposition 2.1 imply that 


o;(£ 4 )= ^ (2.22) 

0esp(E‘) 0esp(E*) 

4»0 


where p^ = w(l 0 (£ 4 )), the probability of measuring the mean entropy production rate to be <f> over the 
time interval [ 0 , 7]. 


2.3.2 Evans-Searles Symmetry 

As stated in Section 1, the first numerical evidence for a fluctuation relation was obtained in [ECM] in 
the context of deterministic classical N particle systems. There, they studied the non-equilibrium steady 
state of a fluid under an external sheer and conjectured that 


Pt{~A) tA 

PM) 


(2.23) 


where Pt (A) is the probability of measuring A for the average dissipation of power over the time interval 
[0, i]. Evans and Searles [ES] were the first to prove that such a relation holds, and we call this relation 
of this sort the Evan-Searles fluctuation relation. 

The direct quantization of the relation (2.23) to our current setting is 

P% = e-“V;. (2.24) 


This equality implies that the Evans-Searles functional 

ES t (a) = log W (e-“ 4S ‘) = log £ e~ at %, (2.25) 

0esp(S t ) 


has the symmetry 


ESt(a) = ES t (l - a). 


(2.26) 


We remark that we sometimes use the notation eES ,t(®) = ESt(a). It is easy to see that (2.26) fails in 
general unless the system is in a steady state. More precisely, we prove the following: 



Proposition 2.2. The functional ESt(a) has the symmetry (2.26) if [H,ui\ = 0. If [H, u>] A 0, then the 
fluctuation relation (2.26) can hold only if t is in a countable set described below. 

Proof. If [H, ui\ = 0 then ESt(a) = 0 and so (2.26) holds. Conversely, suppose the fluctuation relation 
holds for some t. Then 0 = ES,(0) = ESt(l) = logw(e -ts ) and therefore w(e -tl: ) = 1. Furthermore, 

w(e“ tE ‘) = tr (e log u e log _log “). 

By the Golden-Thompson inequality (see Corollary 2.3 of [JOPP]), 

tr(e Iog “e logw - t_Iog “) > tr(e log “‘) = 1, 

and equality holds only if lo and uj-t commute. Let uj = Xk\ek)(ek\, with {e k }k orthonormal. Consider 
the functions 


(2.27) 

(2.28) 


fjk(t) = (ej,e ltH we ItH e k ). (2.29) 

As analytic functions, they either vanish identically or have an isolated, countable set of zeros. Further¬ 
more, oj-t and uj commute iff fjk(t) = 0 for every (j. k) e A := {( j , k ) : A j A Afc}. If there is at least one 
(j, k) e A s.t. fjk(t) is not identically 0, then the set 

{t:f jk (t) = 0\/(j,k)eA} (2.30) 

is countable, and t must be an element of it if the fluctuation relation holds. If all of these functions 
vanish identically then we have that w_t and u> must commute for every t. It follows from the formula 

= i < 2 - 31 > 

\z-E\=e 

which holds for self-adjoint A and e small enough that the spectral projections of u> and commute. 
From the Lemma preceding Theorem XII.8 in [ IS4] the projection-valued functions 

Qjk(t) = 1a, (w)lA fc (w-t) = lA, (u)e itH l Xh (u;)e- itH (2.32) 

have constant kernel (note that it is projection-valued due to the fact that 1^. (cu) and l\ k (w_ t ) commute). 
We conclude that Qj k {t ) = Qj k { 0) = 0 if A, A \ k and that 


IaM = l\(u})e itH ( J] IvH ] e~ HH = lA( W )e i(ff l A (a;)e 

i A'esp(u;) 


-i tH 


Y, lv(w) ] = l A (w-t), 

, A'esp(cj) 


(2.33) 


from which we conclude ui = 0J-t for every t. Differentiating this at t. = 0 yields [H, uj] = 0. 

Therefore, if [H, uj] A 0, then the functions fjk{t) cannot vanish identically for every (j,k) £ A, and 
so t, must be in the countable set described above if the fluctuation relation holds. □ 


Remark 2.1. The proof of the above proposition also yields that ESt(l) > 0 unless either [H, u] = 0 or 
t is in the countable subset described above. 


Proposition 2.2 tells us that the Evans-Searles functional ESt(a) will fail to have the Evans-Searles 
symmetry (2.26) except in trivial cases. In the sequel we will discuss a natural choice of entropic func¬ 
tionals which satisfy the Evans-Searles symmetry. 
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2.3.3 Entropic functionals derived from full counting statistics 

In this section we introduce the concept of the full counting statistics associated to a repeated quantum 
measurement protocol of the entropy flow of our system. The full counting statistics will lead us to define 
an entropic functional that satisfies the Evans-Searles symmetry. 

We define the entropy observable as 

S = - logw = (3 L H L + p R H R + Z, (2.34) 


with Z = tr (e~P LHL -P RHR ). Recall that we have set /3c = 0. We have S t = — logw_ t , and Proposition 
2.1 implies that tYf = St — S. 

According to the postulates of quantum mechanics, a measurement of the entropy of our system at 
t = 0 will yield the eigenvalue s of the observable S with probability u;(l s (S)). After this measurement, 
the system is in the reduced state 


wl s (5) 

w(MS))' 

At a later time t, the system has evolved and is in the state 

itH 

U>± s (S) 


(2.35) 


(2.36) 


A measurement of the entropy at time t will yield the eigenvalue s' of S with probability 

itH 


tr e 


wig (5) 

The joint probability of these two measurements is 

tT(e~ ltH u}± s (S)( 


itH 


tAS)j . 

(2.37) 

MS))- 

(2.38) 


The mean rate of entropy change is <f> = (s' — s) /t, and we have derived its probability distribution, which 
is 

P M) = ^ tr(e _itJ? wl s (S)e itH t s > (S)). (2.39) 

s,s r 

s' —8 = t<f> 

The above discrete probability measure is called the full counting statistics for the entropy change over the 
time interval [0, t\ operationally defined by the specified measurement protocol. The entropic functional 
associated to the full counting statistics is, for a e R, 

FCSt(a) = log£ e~ at *P t (<j>). (2.40) 

4 > 


It is easy to compute, 

FCS t (a) = log]T ^ e- Q ( s '- s )tr(e- itff W l s (^)e itff l s ,(5)) 


s' —S = t(j) 


= log tr e~ a< ' s '~ s ^e~ itH uit s (S)e itH l s '(S')^j = logtr (w}-V) . (2.41) 


With this identity we can verify that the Evans-Searles symmetry holds: 


FCSt(a) = logtr(w t 1_ “w“) = logtr(w i_a w“ t ) = logtr^tW 1-0 ) = logtr^w 1 "") = FCS t (l - a). 


(2.42) 
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The second to last equality follows from an application of the time reversal. This identity implies that 

FCS(a) - FCS(1 - a) = £ e _at * [P t (0) - e^P t (-0)] = 0. (2.43) 

0 


Since this holds for any a e [R, we conclude that 

P t(-0) = e-^P t (^). (2.44) 

Compare this equality with (2.23). This equality says that the probability of measuring a decrease in 
entropy is exponentially suppressed compared to the probability of measuring an increase. 


2.3.4 Entropic pressure functionals 

The definition of FCSi(a) can be generalized. For p > 0 and a e R we define the entropic pressure 
functionals e p ,t(a) by 

I log tr ( to ~p~ OJt to ~~p 

e P ,t(a) = \ V 

[logtr (e( 1_a ) lo 8“+ al °g w ‘) ifp = oo. 

In particular, e 2 ,t(a) = FCSt(a). In the following proposition we summarize the basic properties of 
these entropic functionals. 

Proposition 2.3. The entropic functionals satisfy 

(i) e p , t (0) = e Pjt (l) = 0. 

(ii) The functions R 9 a —* e Pi t(a) are real-analytic and convex. 

(iii) e Pjt (a) = e Pjt (l - a) and e Pit (a) = e Pi _ t (a). 

(iv) The function ]0,oo] £p-> e Pi t(a) is continuous and strictly decreasing for a ^ 0,1, and linip^oo e P: t (a) 

&00,t(&) ■ 

(v) ep t (0) = —tuj(S t ),e' pt ( 1) = tui(Tf). In particular, these derivatives do not depend on p. 

(vi) e 2 , t (a) = e F cs,t(a) and 


if 0 < p < oo, 


(2.45) 


t (0)=ES"(0)= f f w((tr s - w(cr s ))(cr u - w(cr u )))dsdM. 
Jo Jo 


(2.46) 


(vii) eoo ,t(a) = max p S(p,uj) — atplYf), where the max is taken over all states p > 0. 

Remark. We remark that (vii) characterizes 600 , 4 ( 0 ) as the direct quantization of the variational charac¬ 
terization of the Evans-Searles functional in classical non-equilibrium statistical mechanics. 

Proof, (i) These identities follow directly from the definition of e p ,t(a). 

(ii) Real analyticity follows from the strict positivity of to and the analytic functional calculus (i.e., the 
function f(x) = x a is real analytic on (0,oo) and so f(A) is analytic for any strictly positive matrix). 
Corollary 2.4 of [JOPP] implies convexity. 

(iii) The second equality is immediate from an application of the time reversal. Unitary invariance of 
trace norms and the identity ||AB|| = ||BA|| yield 


to/to p 


e~ ltH ujp e ltH to 1 


OL ±1 

ojp e lu p e 


-itH 


UJ p UJ-t p 


= 

1 — a a 

UJ-t P OJP 

p 



(2.47) 
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The first equality then follows from an application of the time reversal. 

(iv). Continuity is clear: for nonnegative matrices, p —* A x / p is a continuous matrix valued function, 
and the trace is a continuous function from matrices to C. Therefore p —* e p t is continuous as it is a 
composition of continuous maps. The remainder of the claim follows from Corollary 2.3 of [JOPP] and 

the fact that we can write e Pjt (a) = log "~ St ■ 


We will prove (v) once we have introduced the Jordan-Wigner transformation (see the remarks fol¬ 
lowing the proof of Lemma 2.10) illustrating the power of the Fermi gas representation. For (vi), direct 
computation yields 


S " t (0) = W ((fE‘) 2 )-(cc (tE*)) 2 

n w(cr s cr u ) - ui(a s )uj(cr u )dsdu 
) 

n w((cr s - U}(a s ))(a u - u(a u )))dsdu. 
) 


(2.48) 


We refer the reader to [JOPP] for the proof of (vii). This property, like the others in this proposition are 
algebraic in nature and hold for general finite dimensional quantum systems. □ 

We close this section with the remark that the variational characterization of e^ ,t(<u) given in (vii) is 
the quantization of the variational characterization of the finite time Evans-Searles functional in classical 
non-equilibrium statistical mechanics. 


2.4 The Jordan-Wigner transformation 

The basic tool in the study of the XY chain is the Jordan-Wigner transformation which dates back to 
[JW]. The Jordan-Wigner transformation maps the XY chain to a free Fermi gas, and associates to the 
Hamiltonian of the XY chain a Jacobi matrix. We assume the reader is familiar with the concept of second 
quantization, and recall here some definitions for notational purposes. For a pedagogical introduction we 
refer the reader to Chapter 6 of [JOPP]. 

Let 1C be a finite dimensional Hilbert space. For n > 1, the n-fold antisymmetric tensor product is 
denoted r n (AC) and is the subspace of the n-fold tensor product of AC with itself, denoted /C® n spanned 
by vectors of the form ipi a ... a ipn^i £ K, where 


ipi a ... a ip n 


—j= Y! si gn(7r)VV(i) ® ® V’Tr(n), 

Vn! 4s n 


(2.49) 


and S n is the symmetric group on n letters. By convention, To(/C) = C. Note that r„(/C) = {0} for 
n > dim AC. For A £ O/c and n > 1, r n (H) and dr„(H) are the elements of Or n (A:) defined by 

r n (A)(^i a ... a 4> n ) = Ail; l a ... a Ail> n 

dr n (H)(^i A ... A ip n ) = Al/; 1 A ... A ip n + ... + Vu A ... A Aij; n . (2.50) 

For n = 0 we set To (4) to be the identity map on To (AC) and set dro(A) = 0. The Fermionic Fock space 
over AC is defined by 


P(/C) 


dim K. 

© r„(/c). 


n =0 


We define for A £ O/c the elements T(4) and dr(H) of Or(iC) 


(2.51) 


r(4) 


dim /C 

© Pn(4), 


n=0 


dr(4) 


dim /C 

© dT n (A). 


71 = 0 


(2.52) 
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We have for A,Be Ok and A e C, 


T(A)* = T(A*), dr(A*) = dr (A)*, 

r(AB) = r(A)T(B), dr(A + a b) = dr(A) + Adr(R), T(e A ) = e dr(A) . ( 2 . 53 ) 

Additionally, one has 


[dT(A),dr(R)] =dT([A,B]), (2.54) 

and for A invertible r(A) 1 = r(A _1 ) and 

T(A)dT(B)r(A~ 1 ) = dT(ABA~ 1 ). (2.55) 

Let fl = 1 e r 0 (/C). The annihilation and creation operators for if) e /C are elements of Or(iC) and are 
denoted by a(ip) and a* (ip ), respectively, and are defined by 

a*(lp)Ll = Ip, a* (ip) ('ll) l a ... A lj) n ) = ij) A Ipi a ... a ip n (2.56) 


and 


a(ip)Ll = 0, a(ip)ip i = 

n 

a(ll))(ll> 1 A ... A l/) n ) = ^](-l ) 1+J (l/),l/)j')Ki’l A ... Al//f A ... A lj) n . (2-57) 

j = 1 

The maps ip —» a*(ip) and ip —» a(V’) are respectively linear and antilinear, and a(ip)* = a*(ip). They 
obey the canonical anticommutation relations (CAR) 

{a(ip),a((p)} = {a*(ip), a*(<p)} = 0, {a(ip), a*(<p)} = (ip,(p)>K lr(AC)- ( 2 -58) 

for any ip and <p in /C. Here {A, B } = AB + BA denotes the anticommutator of two observables. We will 
use the notation a& to refer to both a and a * simultaneously. 

Given finite dimensional Hilbert spaces K. and 1Z, a representation of the CAR over K. on 1Z is a pair 
of maps from K, to On denoted by 

iP^b(iP), iP->b*(iP), (2.59) 

where the first is antilinear and the second is linear, b(ip)* = b*(ip) and (2.58) is satisfied with a# replaced 
by b&. The representation is called irreducible if the commutant of the set {b#(ip),ip e Kf) is trivial, that 
is, if 

{Be On - [.B , b*(ip)] = 0 \/ipeK,} = Cl n- (2.60) 


For finite dimensional K, and 1Z we have the following characterization of irreducible representations of 
the CAR over K. on 1Z (see Exercise 6.2 of [JOPP] and also [L]). 

Proposition 2.4. Let K. be a finite dimensional Hilbert space and ip — » b#(ip) be an irreducible represen¬ 
tation of the CAR over /C on the finite dimensional Hilbert space H. Then, there exists a unitary operator 
U : T(/C) — » 1Z such that Ua&(ip)U* = b^(ip) for all ip e 1C. Moreover, U is unique up to a phase. 

Consider an XY chain confined to a finite interval A c Z, as defined in Section 2. The Jordan-Wigner 
representation is an irreducible representation of the CAR over hv = ^ 2 (A) on 77 a- The resulting unitary 
operator Fjw : 77a —» T(flA) guaranteed by Proposition 2.4 is called the Jordan-Wigner transformation. 
In what follows we construct the Jordan-Wigner representation (see also [JOPP], although the model 
there is slightly less general). 

Consider the spin raising and lowering operators 


a 


(+) = 



r(~) 



(2.61) 
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We have that cri ±) = (<ri 1 ' ) ± icri 2 ^)/2. Moreover, a x ^ and o4 +) = a x obey 

Wi +) > <4 +) } = {°i _) > cr i _) } = 0 {°i +) > <4 _) I = ■ 


(2.62) 


If A = {a;} where just a single lattice site, then the maps £ 2 ({x}) 3 a —* aa e Ou {x) and £ 2 ({x}) 3 

a —» e Ou {x) would indeed define a representation of the CAR over £ 2 {{x}). However if A contains 

at least two distinct sites x ¥= y, then this does not directly generalize, as the raising and lowering 
operators at different sites will commute, and not anti-commute. For example, if A = [0,1], then the 
maps £ 2 (A) 3 (ao,ai) —» aoOo + ^ + e Oa and £ 2 (A) 3 (ao,au) —» ao&o ^ + aqcrj - '* e Oa do not 

define a representation of the CAR. 

Let A = [A, B\. A computation using the fact that ( cr £^) 2 = 1 yields, for T x , S v e Oa, 


{{v£\T x }(T { *l l ...af± l Sy x<y, 
{va ] ■■■°x-i t x,Va ) -"vf-iSy} = j {T x ,Sy} X = y, 

l{4 3) > Sy^h-VxliTx x>y. 

Define for x e A\{A}, 

b — fT*' 3 '* (T^ b* — (T*- 3 - 1 <7^ /y(S~) 

o x — &A ■" a x-l <J x > ° x — a A 


(2.63) 


(2.64) 


and bA = <? ^ \ b* A = cr^ + \ We have 

Proposition 2.5. The maps £ 2 {A) 3 a—* b(a) := Yj X a xb x e Oa and £ 2 {A) 3 a —» b*(a) := Yj x a xb x e 
Oa form an irreducible represenation of the CAR over £ 2 ( A) on TLa- 

Proof. That b # obey (2.58) follows immediately from (2.62) and (2.63). Clearly they are linear (resp., 
antilinear) and 6(a)* = b*(a). To see that the representation is irreducible, we require the following 
lemma, which is Proposition 6.4 in [JOPP]. 

Lemma 2.6. A representation if —> b^(i/j) of the CAR over the finite dimensional Hilbert space 1C on TZ 
is irreducible iff the smallest *-subalgebra of On containing the set {b^(if) |t/> e JC] is On¬ 


Now , define 


14 = 


1 

n ye [A ,4 2b * b * - !) 


if x = A, 
else. 


We claim that the following holds: 


4 1) = v x (b x + b* x ), 4 2 > = iv x (b x -b*), 


4 3) = 2b x b* - 1 . 


(2.65) 


( 2 . 66 ) 


Since the Pauli matrices form a basis for C 2 and all polynomials in the operators b # are in the smallest 
*-subalgebra containing {&#(?/>) \f) e £ 2 (A)}, the relations (2.66) show that this *-subalgebra is indeed all 
of Oa- 

Finally, to prove (2.66) note that 2b*b y — t = 2cr]) + Vy _ ’ > — t = cry 3 ' 1 . The other two relations follow from 

this and the equations (ai 3 '*) 2 = 1 and a x ~^ = ( ai 1 ' 1 ± \cr^)/2. This completes the proof of Proposition 
2.5. □ 

As previously noted, this proves: 

Proposition 2.7. There exists a unitary operator Ujw : TLa r(f)A), called the Jordan- Wigner trans¬ 
formation satisfying 

U Jw a^U^ = S x (a x + a* x ), = iS x (a x - a* x ), t/jw<4 3) t/jw = 2a* x a x - 1 (2.67) 


where 


t x = A, 

Yl y e[A,x[^ a x a x ~ 1 ) e ^ se - 


( 2 . 68 ) 
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Consider now the open XY chain confined to the interval A = as defined in Section 2.2, 

with left part A; := [— M, — N — 1], central part A c := [— N, AT] and right part A r := [N + 1 ,M], The 
Jacobi matrix h associated to our XY chain is a bounded operator on ^ 2 (Z) given by 

(. hu)(n ) = J n u(n + 1) + J n -\u(n — 1) + v n u(n), (2.69) 

for u e £ 2 (Z). For any subset A of the integers we denote for k e A the element 5 k e £ 2 (A) which is 1 at 
k and 0 elsewhere. 

Recall that 1 )a = £ 2 (A) and we define \)i/ c / r = £ 2 {Ai/ c / r ). They are the single particle Hilbert spaces 
for the free Fermi gas on £ 2 {A/ A l / c / r ). Let h\ and hi/ c / r be the restrictions of the Jacobi matrix h to 1 )a 
and fyi/ c / r , respectively. The h\ and h[/ c / r are the single particle Hamiltonians. For simplicity of notation, 
we identify the operator hi acting on hi with the operator hi © 0 © 0 acting on 6 = (p © t) c © hr, and we 
make similar identifications for h c and h r . Note that 

h\ = ho + v, (2.70) 

where ho = hi + h c + h r and v = v r + vi with 

vi = J-n-i(\5-n-i)(5-n\ + |5_j\rX<5_jv_i|), (2-71) 

v r = Jn{\6n+i)(8n\ + | 5n)(6n+i\)- (2.72) 

Up to an irrelevant additive constant, we have the following identities which are an immediate conse¬ 
quence of Proposition 2.7: 

UjwHi/c/rU^ = dT(hi /c/r ), UjwVi/ r U^ = dr(u ; / r ), UjwHU^ = dr(/i). (2.73) 

For the fluxes we have 

UjW^lUjw = —iJ-N-lJ-N-2(a-N a -N-2 ~ U-N- 2°-Jv) _ iJ-N-1 A-JV-l («* N a -N-1 ~ O* jv-1°-Jv) 

(2.74) 

CJjwd’rCJjw = — iJNJN+i(dN a N+2 — a% +2 ON) ~ iJ jvAjv+i( a^ajv+i — a^, +1 aAr). (2-75) 

In the fermionic representation, it is easy to see that the XY chain is time reversal invariant. If j is 
the standard complex conjugation on £ 2 (A), then conjugation with the anti-unitary operator 

0 = U^T(j)U JW (2.76) 

is a time reversal for the confined XY chain under which the initial state is invariant. 

For the remainder of the paper we will work only in the Fermionic representation of the XY chain. 
We will slightly abuse notation and write Ha = dr(/iA), H l / c / r = dT(hi/ c / r ). Note that the proofs of 
Propositions 2.1, 2.2 and 2.3 carry over in this representation without change. 

2.5 Basic Formulas 

In this section A = [— M, M] is fixed and we will drop the respective subscript and write h for h\, H for 
H\, f) for [)a, etc. We would like to derive some basic formulas for the entropic functionals. First, let us 
record a few basic identities that will prove useful in this section and the sequel; for the proofs, we refer 
the reader to [JOPP]. We have, 


tr(r(A)) = det(l + A), 


(2.77) 


which holds for any linear A : fj —» fj. 

Some of the basic algebraic properties of the annihilation and creation operators are summarized in 
(see Proposition 6.2 of [JOPP]), 

Proposition 2.8. The following holds: 
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(2.78) 


(i) ||a # ('0)|| = HV’II for any if, 

(ii) For any A e 0^, 

r (A)a*('if) = a*{Axj))T{A), T {A*)a(Aip) = a(t/j)r(A*), 


(iii) If U is unitary, 


(iv) For any A e O r,, 


T{U)a*{i/j)T{U*) = a*{Uif), 


[dr^aXVO] = a*{Aif), [dr(^),a(V’)] = -a(A*ij>), 


(2.79) 


(2.80) 


(v) a*{i/j)a{(/)) = dr(|^><(/>|). 

Define k = —pihi — /3 r h r and fc t = e lth ke~ lth . The initial state of the XY chain is mapped to the 
density matrix, which we also denote by uj, 


,-faIh-prHr 


D dr(e fe ) 


UJ = 


tr(e P lHl P rHr ) det(l + e k ) 
The state uj is quasi-free with density T = (1 + e -fc ) -1 . That is, 

1 / T 


UJ = —r 


Zj* V 1 — T 


where 


Zt = tr T 


T 


t-T 


We have the following (see Proposition 6.6 in [JOPP]): 

Proposition 2.9. Since uj is quasi-free, 

(i) If ,ip n e t), then, 

u(a*((f n )...a*...afifm)) = S nm det^i, T<j>j)\, 

(ii) w r (r(A)) = det(l + T(A - 1)), 

(iii) W7’(dr(yl)) = tr (TA). 

By the definition of the entropic functionals and the identity (2.77) we have for p < oo, 

det(l + ( e ( 1_ “) A: / p e 2 “ fc - t ' /p e^ 1_a ^/ p ) p / 2 ) 


e P ,t(a) = log- 


We also have 


and, 


eoo ,t(a) = log 


det(l + e fc ) 
det(l + e (i-«)fe+«fe-t) 


ES t (o) = log 


det(l + e k ) 

det(l + e fe / 2 e «( fe t- fc ) e fc / 2 ) 
det(l + e k ) 


(2.81) 


(2.82) 


(2.83) 


(2.84) 


(2.85) 


( 2 . 86 ) 


(2.87) 
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Define for p e]0, oo[, 


/Cp, t (a,«) = l+(e( 1 

-o)k tu /p e 2o:fe_ t(1 _ u ) /p e (l-a)k tu /p\ P ^ 

e (l-a)k tu /p + h c ( 

V - 

- / 

(2.88) 

where h.c. stands for the hermitian conjugate of the first term. Define 


A^oD,i(cfcj u) 

— g (1 ot)kty otk_ t(i_ u ) ^ 

(2.89) 

and 



/Ces ,t(o:,u) = 

_ ^ _j_ e -<x(k t(1 _ u )-k-tu) e -k- tu \ 

(2.90) 

We have, 



Lemma 2.10. For p e] 0, oo], 

rot rl 


£'P,t(p8) t 

d 7 dwtr(/C p> t( 7 ,w)i[fc,ft]), 

(2.91) 


'o Jo 


and 

-*CA 


ESt(a) = t 

d7 dutr(/C E s,*(7>u)i[fc,/i]). 
o Jo 

(2.92) 


We will need the following result of [HP], of which we provide a proof: 

Lemma 2.11. Let f(z) be an analytic function in a disc of radius r centered at some A e R. Let F{pt) 
be a Hermitian matrix valued function in a e (—<5, 6) for some 8 > 0. If F(a) is differentiable at a = 0 
and the eigenvalues of F{ 0) are contained in (A — r, A + r), then 


= tr(/'(F(0))F'(0)). (2.93) 

ot=0 

Proof. By continuity of F(ct) at a = 0 and min-max principles, the eigenvalues of F(ct) are all contained 
in (A — r, A + r) for |ck| < <5o where <5o > 0 is a constant possibly smaller then 8. By the analyticity of f(z), 
the function a —> /(F(a)) is defined by a norm convergent power series f(F(a)) = Ym=o c n(^{ a ) ~ A) n 
for | a | < Sq. By this fact and the continuity of the trace, 


— tr/(F(a)) 

da 


tr f(F(a)) - tr/(F(0)) 


^E 

n=C 

00 

= E 

n=C 

00 

^ E 


-tr (/'(F(0))F'(0)) 


tr ( (m-w-mo) + xr _ n{m _ A) »-i no) 


tr 


V a 

(F{a)-F{ 0) 


\ a 


P n {a) - n(F(0) - A) n ~ 1 F'(0) 


n =0 


tr P n (a) 


F(a) - F(0) 


\ 00 

P'(0) + E M l trF '(°) ( P 4«) - P «(°))l> ( 2 - 94 ) 

' n =0 


with P„(a) = (F(a) - A) 71 " 1 + (F(a) - A) n ~ 2 (F(0 ) - A) + ... + (F(0) - A) ra_1 for n > 1 and P 0 (a) = 1 
(here we have used cyclicity of the trace to obtain the equality in the above computation). Continuity 
of F gives us the estimate that, by taking 5o smaller if necessary, there is an ao e (—(So,<Ao) so that 
|[P n (a)|| < n H.F(a'o) — A|[" holds for every |or| < 8q/2. By differentiation, we conclude that for any 
e > 0 for we have |[P n (a) — P„(0)|| < en(n — 1) ||F(ao) — A||" -2 for every |a| < Sq/2 (again taking <J 0 
smaller if necessary). 

Since Jf ln n(n — l)\c n \ ||F n (ao) — AH" - converges, these estimates together with the differentiability 
of F at a = 0 show that the last line of (2.94) goes to 0 as a —» 0. □ 

This lemma implies, 
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Corollary 2.12. Let F(a ) be a differentiable function for a e IR taking values in Hermitian strictly 
positive matrices on f)- Then for any p e]0,oo[, 


— trlog(l + F(a) p ) = ptr ^(l + F(a) p ) 1 F(a) 1 F'(a)^j . (2.95) 

Let F(a) be a differentiable function for a e IR taking values in Hermitian matrices on f). Then, 

trlog (l + e F(a) ) = tr ((1 + '(a)) . (2.96) 

We may now proceed with the computation of the basic formulas for the entropic functionals. 

Proof of Lemma 2.10. We will complete the computation for 0 < p < oo. Recall that for matrices 
log det A = tr log A. By this fact, 

e p , t (a) = trlog + ^e (1 - a)fe / p e 2Q/c -‘/ p e (1 - a)fc / p ) P/2 ^ - trlog(l + e k ) 

= J d 7 -^ trlog fl + ( e (1 - 7)fe/p e 27fc -* /p e (1 - 7)A: / p ) P/2 ^ . (2.97) 

Let Ffy) = e( 1_7 ) fe / p e 27fe_t / p e( 1_7 J fc/,p . By (2.95), the derivative in the above integral equals 


tr ((1 + F( 7 )- p / 2 )- 1 F( 7 )- 1 (e (1 - 7)fc/p (fc_ t - k)e 2 ^ k - tlp e {1 ~^ klp + e {1 ~^ klp e 2lk - tlp fk_ t - k)e {1 ~^ k/p jj 
= tr ((1 + F[^)- p / 2 )~ l (e^-^Pfk-t - k)e~ (l ~ l)k/p + e" (1 ~ 7)fc / p (fc_ t - k)e^ Wp ^j ) , (2.98) 

where the last equality follows from cyclicity of the trace. We compute next, 

k_ t — k = f ds-^-k_ s = i f d se~ ish [k,h]e ish = it f dse~ isth [k, h]e isth . (2.99) 

Jo ds Jo Jo 

Inserting (2.98) and (2.99) into (2.97) and some simple algebra yields the claimed formula. 

The computation for p = go follows the same strategy. Taking F(a) = (1 — a)k + ak-t in (2.96) 
yields, 


-oo,tt 


*) = J d 7 tr (( 


1 + e 




)- 1 (fc_ t -fc)) 


= t 

= t 


m 

M 


du tr 


du tr 


1 _|_ g—ig—i«^j 
1 + e -(l-7)fct u -7fe-t(l-„))-li 


[k,h]e iuth ^j 

[k, ft]) 


( 2 . 100 ) 


which is the formula in question. The computation for ESt(a) is similar and can be found in [ ]. □ 

Evaluating the derivatives computed in the proof of Lemma 2.10 at a = 0 and using part (iii) of 
Proposition 2.9 to evaluate 5 , (wt|o;) in terms of k and k_ t yields (v) of Proposition 2.3. 


3 The extended XY chain: thermodynamic limit 

Recall that our open XY chain restricted to A = [— M,M ] consists of a left part A; = [— M, —N — 1], a 
center part A c = [—AT, AT] and a right part A r = [AT + 1 ,M], We now consider the limit \—M,M] —* Z 
in which we keep the central system A c = [— N, AT] fixed and take M —* oo. Recall that we are working 
exclusively in the Fermionic representation of the XY chain. To denote the dependence of the various 
objects under consideration on the size of A we use the subscript M and write, for example, Om, 1) m, hM, 
etc., whereas before we wrote Oa, F)a, Aaj etc. 


18 


The algebra of observables of the extended XY chain is denoted by O and is the norm closure of the 
local observables 


Oioc = (J O m , (3.1) 

M 

where we have identified Omi with the appropriate subalgebra of Om 2 for Mi < M 2 . Recall here that 
O m = and that 

T(fiM 2 ) = r (^ 2 (Am 2 \AmJ ©^ 2 (AmJ) = r (^ 2 (Am 2 \Ami)) <8>r (t? 2 (A Ml )) ■ (3.2) 

For the dynamics of the XY chain we have, 

Proposition 3.1. For any A e O\ oc the limit 

tHA) = lim Tm(A) = lim e ltHM Ae~ itHM (3.3) 

M— >oo M— >oo 

exists in norm and the convergence is uniform for t’s in bounded sets. Furthermore, r* uniquely extends to 
a strongly continuous group of *-automorphisms of O. This group is called the dynamics of the extended 
XY chain. 

In the above proposition, strong continuity means that the function t —*■ t*(A) is norm continuous for 
every As O. 

For the inital state of the XY chain, we have, 

Proposition 3.2. Let ojm be the initial state in (2.81). Then for A G 0\ ac the limit 

u)(A) = lim ujm(A) (3-4) 

M— >oo 

exists and ui uniquely extends to a state on O. This state is called the initial state of the extended XY 
chain. 

The extended XY chain is described by the C* dynamical system (0 ,t*,u;). This system is time 
reversal invariant, as the time reversal defined in Section 2.4 which is initially defined for every A e 0\ oc 
extends to all of O by the bounded linear transformation (BLT) theorem [RSI]. We denote this time 
reversal by ©. 

Proof of Proposition 3.1. We first prove that the norm convergence of the limit is uniform for bounded t. If 
A G Oi oc then A is a finite sum of polynomials in the annihilation and creation operators {a^(if)\if g t Jm 1 } 
for some M'. It therefore suffices to prove the uniform convergence when A = a* (cf n ).. .a* ((fi)a(ifi).. .a(if m ) 
For M > M', we have by Proposition 2.8(iii), 

e itdr(h M ) a *^ .. a (il, m ) e - itdr (hM) = a * {e lthM (f n )...a* (e lthM 4>\)a(e lthM , f)x)...a(e lthM xfrn). 

(3.5) 

Since \\a#(e lthM ip)\\ = ||^|| for any if, it suffices to prove that, for any if G f )m’, the norms of the sequence 
e ithM^p are urL iformly Cauchy for bounded t. However this follows directly from the fact that Lm —» h 
strongly in £ 2 (Z) (recall h is the full Jacobi matrix defined in (2.69)) and the estimate, 

n -* oo 

| \e ithM if - e ith if | < ^ -\t\ k | (h k M - h k )if || + \\if\\ - \\h\\ k + \\h M \\ k • (3.6) 

k =0 ' k=n +1 

Since the ||/im||’s are uniformly bounded, the second sum can be made smaller than any e > 0 for n large 
enough, uniformly for bounded t, and then since h^ —» h k strongly for every k , the first sum is also 
under control. 

Since ||e ltffM || < |]A|| holds for every M and A G 0\ oc , t 4 extends uniquely to all of O by the 

BLT theorem [RSI]. Since satisfies the group property and the ^-automorphism properties on O i oc , 
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it satisfies the same properties on O. Strong continuity on O\ oc follows from an e/3 argument using the 
fact that the convergence is uniform for Vs in bounded sets. Strong continuity on all of O then follows 
immediately from a second e/3 argument. □ 

Proof of Proposition 3.2. Let A e O i oc be given. As in the proof of Proposition 3.1, it suffices to consider 
A = a*(<f„)...a*(<j>i)a(if>i)...a(if>m), with the (ffs and ipj’s in \)m’ for some M'. By Proposition 2.9(i) we 
have for M > M ', 


uj T {a* (4> n )...a* ((fi)a{'il) 1 )...a('>p rn )) = 5 nm det[<^i, T M ^j)], (3.7) 

with Tm = (1 + e~^ lhl ’ M ~^ rhr ’ M )~ 1 . The existence of the limit (3.4) is then immediate from the strong 
convergence of hi/ r M to hi/ r . The estimate \u>m{A)\ < ||A|| which holds for every M implies that u> 
extends uniquely to all of O. □ 

The observables <&i/ r and a are in 0\ oc . Let 

1 r* 

= - cr s ds, (3.8) 

t Jo 

and 

ES t (a) = logo; (e" atE ‘) . (3.9) 

We prove the following two basic propositions about the observables and entropic functionals in the 
thermodynamic limit. 

Proposition 3.3. We have, 

(i) 

lim Yf M = E* (3.10) 

M—► oo 

in norm and the convergence is uniform for t’s in bounded sets. 

(ii) liniA/_»oo S(u ttM \uM) = ~tw (E*) 

(iii) r* (E -t ) = E‘ 

(iv) E ( = —t* (0(E 4 )). In particular, sp(E ( ) is symmetric w.r.t. to the origin. 

(v) 


lim ESt m{oI) = ESt(a) (3-11) 

M —>oo 

and the convergence is uniform for t’s and a’s in bounded sets. 

Proof, (i) The function r^(a) converges to t s (ct) uniformly for s e [0, t] due to Proposition 3.1 and the 
claim follows. 

(ii) The definition of w and (i) imply that —tuiiT, 1 ) = limjv->.oo liniM->oo E^) = liniM->oo — Iwm/E^). 

The second equality follows from the uniform norm continuity of the state u>. By Proposition 2.1, 
— tuiMO&M) = >S'( w t,M|wM) and the claim follows. 

(iii) From Proposition 2.1 and (i) we have E* = liniM->oo = hm.M-.oo t m (^m)- O n the other hand, 
r* (E _t ) = limjv->oo hmM^oo t~m {^~n) an d therefore (iii) follows from the inequality |r^ f (E^ t — E/^)|| < 

II E“* — E-;|, which implies 

lim lim t^(E-*) = lim r^(E^). (3.12) 

iV —>00 M—>oo M—>oo 
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(iv) We have, 


—-r(©(£*)) = lim -^(0(2^)) = lim -t‘(0m(SV)) = lim lim -r^(0 M (^)). (3.13) 

M->oo M->oo M—>oo AT-»oo 

The RHS equals the limit limM^oo -4 f (0M(Sji,;)) which equals limM->oo Em by Proposition 2.1, which 
proves (iv). 

(v) Note that the RHS of (3.11) equals limAr^oo limM—oo log WM(e"“ ffi “) by definition. By the inequality 


Iw M (e _atEtM - e _ “ ts ‘ N )| < 
where the RHS goes to 0 as N, M —* oo, 


^—atY, M g—atE N 


(3.14) 


lim lim logWM(e atEw ) = lim logWM(e 
N ->oo Af—>oo M—>oo 


-atE, 


(3.15) 


This limit is the LHS of (3.11). The desired uniform convergence follows easily by expanding the ex¬ 
ponential on the RHS of (3.15) and the uniform convergence of the (recall the argument using the 
estimate (3.6) in the proof of Proposition 3.1). □ 

Proposition 3.4. For the entropic functionals, we have 

(i) For all a e IR and p e]0, oo] the limits 


e P ,t(a 0 = lim e Pjt ,M(ac), 
M—> oo 


(3.16) 


exist and are finite. 

(ii) e P) t(0) = e Pit (l) = 0. 

(iii) The functions a —» e Pj t(a) are real-analytic and convex in a and jointly continuous in ( p,t,a ). The 
function a —» ESt(a) is real-analytic in a. 

(iv) e p>t (a) = e p - t (a). 

(v) e Pjt (a) = e p>t (l - a). 

(vi) The function ]0, oo] 3 p —» e Pjt (a) is continuous and decreasing. 

(vh) e ; >t (0) = ESJ(O) = - e ; it (i) = -M s*). 

(viii) 

e2,t(0)" = ES"(0) = f f w((cr s - w(cr s ))(cr„ - w(cr„)))dsdM. (3-17) 

Jo Jo 

(ix) As M —* oo, the sequence Pt,M converges weakly towards a Borel probability measure P* on R and, 

e 2 ,t(a) = log f e _to0 d4>t(</>). (3.18) 

Jr 

All the moments of Pt,M converge to corresponding moments of Pt. The measure Pt is the full 
counting statistics of the extended XY chain. 

Proof. Recall the formulas (2.91) and (2.92) and that Iim —>■ h strongly. The also have 

strong limits which we denote /C# jt ( 7 , it), where ff stands for either p e]0, oo] or ES. Note that their 
strong limits are given by (2.88), (2.89) and (2.90), but with the finite dimensional objects replaced with 
their infinite dimensional counterparts (i.e., the kM and hM replaced by k and h which are operators on 
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£ 2 (Z)). Note that kt is given by e lth ke~ lth , and not the dynamics of Proposition 3.1 (the r 4 act on O 
which is a different object than the bounded operators on l 2 (l) of which h is an element). 

Since the commutators [/cm, ^-m] = [fc, h\ do not depend on M and are finite rank, the traces in (2.91) 
and (2.92) are the sum of only finitely many terms. Explicitly, one has (recall that we previously defined 
&ES, t,M( a ) = ES tiM -(a)), 

e#,t,M(a) = t E I dy du(6k, lCft't,Mil, u)i[k, h]Sk) (3.19) 

keA 

where A = [—N — 2, —N + 1] u [N — 1 , N + 2] c Z. In particular, the size of A is finite and does not 
depend on M. It is easy to see that there is a constant C depending only on a and # s.t. 


||£# 1 t,M(7,w)ll ^ C. 


(3.20) 


for all ( 7 ,u) G [0,a] x [0,1] (here, if a < 0, [0,a] := [a,0]). It follows by dominated convergence that, 


lim e #>tjM (a) = t V dy 

M ^°° Jo 


d u lim <4,/C #itj M(7>w)i[ fc i^]4) 
M—> 00 


= f d'l f dit<4,/C #>t (y,M)i[fc, h\5k) 
keA J 0 Jo 

= t\ dy dutr (/C #jt (y,u)i[fc, h]), 

Jo Jo 


(3.21) 


which proves (i). 

The /C Pi t(y, u) are jointly norm continuous in (p, f,u, y) and so the functions (4>^#,t(7> u )4) are 
jointly continuous in (p, t,u, y). Since the sum in the second last line of (3.21) is finite, the functions 
e Pi t(a) are jointly continuous in (p, t , a). Let us be careful in the proof of the analyticity of the entropic 
functionals. We claim that the analyticity of the functions 


a -> <4, /C#,t(7> «)i[fc, h]6 k ) (3.22) 

follows directly from Lemma A.l. From this, the analyticity of the e#y(a) is immediate. If # = 00, then 

£oo, t (a,s) = /(Ei(a, S ,f)), (3.23) 

with f(z ) = (1 + z) _1 and Fi(a,s,t) = e ( - 1 ~ a ^ kta ~ ak - t( - 1 - a '>. It is clear that Fi(a,s,t) > e - I 1+2 l a l)ll fc ll. 
The series expansion for F 3 (a, s, t) converges for every a, s and t, and the coefficients of the power series 
expanded about some ao are dominated by the coefficients appearing infront of the \a — op I in the series 
expansion of e ( 1 + 2 !“-“ol+2|ao|)l|fc||, anc j so Lemma A.l applies. 

If # = P, with p g] 0 , 00 [, then 

>C p ,t{a,s) = F 2 (a, s,t) f(Fi(a, s,t))F 3 (a, s,t) + F 3 (a,s,t)f(F 1 (a,s,t))F 2 (a,s,t), (3.24) 

with f(z) = (1 + z~ p Z 2 )" 1 , F 2 (a,s,t ) = , F 3 = F^ 1 and, 

Ft(a, s, t ) = e {l - a)ktalp e 2ak -^-^ lp e {1 - a)kta / p > e 1 ( 1+4 l“l)ll fc I/p. ( 3 . 25 ) 


The dominating series are obtained by the expansions of e^ 1 + l a a °\+\ a o\) K /p f or p 2 and p 3 , and 
e (i+4|a-ao|+4|ao|)||fe||/p for p^ and so L emma A.l applies. 

For # = ES, we have 


4es, t(a, s) = -F 2 (a,s,t)f(Fi(a,s,t))F 3 (a,s,t) 
with f(z) = (1 + z) _1 , F 2 (a,s,t) = e _fc -* s / 2 , F 3 = F 2 _1 and, 

Fi(a,s,t) = e~ k - tsl2 e~ a ^ 1 - a 'i~ k - ts) e~ k - tal2 > e —(1+2|ct|)l|fc|1 . 


(3.26) 


(3.27) 
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The dominating series are given by the constant e K for F 2 and F 3 , and e( 1 + 2 l“ - “°l + 2 l a °l^ fc ll for F±. This 
completes the proof of analyticity. 

The remaining parts of (ii)-(vi) follow from the corresponding properties of the e# t t,M(&) collected in 
Proposition 2.3. 

For (vii) and (viii) we need the following (see Appendix B in [JOPP]), 

Theorem 3.5 (Vitali’s convergence theorem for analytic functions). Let D( 0, e) c C be the disc of radius 
e centered at 0, and let F n : D(0,e) —♦ C be a sequence of analytic functions such that 

sup |Fk(")| < go. (3.28) 

zeD( 0,e) 
n 

Suppose that the limit 


lim F n (z) = F(z ), 

n—> go 


(3.29) 


exists for every z £ D{ 0, e) n IR. Then the limit (3.29) exists for all z e D{ 0, e) and is an analytic function 
on D( 0,e). Moreover, as n —* go, the derivatives of F n converge uniformly on compact subsets of D( 0,e) 
to the corresponding dervatives of F. 

The arguments that gave the analyticity of the e# t t(a) also apply to the finite volume by 

replacing the infinite volume versions of the £#^( 7 , u) with the finite volume 7 ; u) (the argument 

needs only to be modified by replacing ||fc|| with sup M \\kM\\ < 00 in the bounds appearing there, and 
also by fixing t and allowing the Fj(a, s,-)’s to be indexed by M in their third argument instead of t). 
It follows from Lemma A.l that for each ff there is an ro > 0, which does not depend on M s.t. the 
functions all have analytic extensions to the disc Z?(0,ro), and so by the Vitali convergence 

theorem, Proposition 2.3 and (i), we conclude (vii) and also the equality 


e 2 ,t( 0 )" = ES"(0) = lim cv M ((cr s M - UmWs,m)){v u ,m - WM^M^dsdu (3.30) 

M^coJ 0 J 0 

Similar arguments as in the proof of Proposition 3.3 allow us to conclude that the integrand converges 
uniformly to ui((cr s — Lv(cr s ))(a u — u>(cr u ))) for (s,u) e [ 0 ,<] x [ 0 ,f] which proves (viii). 

To prove (ix), consider the sequence of functions 


J e" Qt *dP t , M (0) = e 


_ p e2,t,M (ot) _ 


= det (l + (1 + e~ kM )~ 1 ( 


det(l + a )^ M e~ l thM ^th M ^ 
det(l + e kM ) 


1 )) 


(3.31) 


which are all entire analytic. The bound | det(l + A)| < where HA^ denotes the trace norm of A, 

together with 


g—a/cMg— ithM^otkM ^ithM _ ^ _ g—afcjUg—is/iMUgafc*f ^Je is/lM ds 

Jo 


(3.32) 


(with v finite rank) imply that for any bounded set BcC, 

sup |e e2 '*’ M(a ) | < 00 . (3.33) 

aeB 

M 


Here we have again used the uniform boundedness of the ||/im||’s. By the Vitali convergence theorem, 
the sequence of characteristic functions of the measures P ty M converges pointwise to an entire analytic 
function, and the convergence is uniform on bounded sets. The existence of and the convergence to the 
weak limit P t follows from the lemma following this proof, and the convergence of the moments follows 
from Vitali’s theorem. This proves (ix). □ 

The following lemma is Corollary 1 to Theorem 26.3 in [Bi] 

Lemma 3.6. Suppose that p n are probability measures with characteristic functions (f> n and that 
limn -,00 4> n {t) = g(t) for each t where the limit function g is continuous at 0. Then g is the characteristic 
function of a measure p which is the weak limit of the p n . 
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4 The large time limit 

4.1 Operator theory preliminaries 

We first collect a few facts about spectral and scattering theory which we will use without proof. The 
reader is referred to [Ja] (especially Section 4.9) for proofs and a more complete exposition. Let A be 
a self-adjoint operator on a Hilbert space /C. The spectrum of A is denoted sp(A) and its absolutely 
continuous part sp ac (4). The projection onto the absolutely continuous subspace of A is denoted l ac (4l). 
For any ip i, '02 £ /C, the boundary values 

(ipi, (A- E± iO)“V 2 > := lim<V>i, (A - E ± ier) -1 g/> 2 ), (4.1) 

exist and are finite for Lebesgue a.e. E £ IR. Whenever we write {ipi, (A — E + iO) -1 ^) we assume that 
the limit exists and is finite. If is the spectral measure of A for ip, then the Radon-Nikodyn derivative 
of its absolutely continuous part is 

di'ty.ac (E) = -{ip, (A- E - iO)~ 1 ip')dE. (4.2) 

7r 

Let Si := S-n -i and 5 r := 5 n+ i- We denote by vi/ r the spectral measure of hi/ r for 8i/ r . By (4.2), 

4^/ r ,ac = -Fi/ r {E)dE (4.3) 

7 r 

where 

F l/r (E) = Im G l/r (E), G l/r (E) = (5 l/r , (h l/r -E- iO)" 1 ^) (4.4) 

We will also denote, for z £ €\IR, the resolvent Gi/ r (z) = {Si/ r , ( hi/ r — z)~ 1 Si/ r }. 

By the spectral theorem, we may identify f)ac(^o) with L 2 ([R, d^; iac ) © L 2 (IR, dry iac ) and ho (fj ao (ft. 0 ) 
with the operator of multiplication by the variable E £ IR. The set S;/ rac = {E £ IR| Fi/ r (E) > 0} is 
an essential support of the absolutely continuous spectrum of hi/ r . Furthermore, let J-n- i = Ji and 
Jn = J r - We set 


£ — Sz, ac n Yj r ac . (4.5) 

Let us recall basic facts from trace-class scattering theory, the proofs of which can be found in [RS3] , 
Section XI.3 (see especially Theorem XI.8). The basic existence result is, 

Theorem 4.1 (Kato-Rosenblum theorem). If A and B are bounded self-adjoint operators and A — B is 
trace-class, then the wave operators 

w+(A,B) = s-lirn e itA e~ itB l ac (R) (4.6) 

£—► + 00 

exist and are complete. That is, 

Ranui_(4,R) = Ranizj + (4., B) = t ac (A). (4-7) 

Moreover, sp ac (4) = sp ac (I3) and the essential supports of the absolutely continuous spectrum of A and 
B coincide. 

Note that If the wave operators exist and are complete, then the scattering matrix 

s{A,B) = w+{A,B)*w-{A,B), (4.8) 

is a unitary operator h ac (B). We have also, as a consequence of completeness (see Propositions 1 and 3 
of Section XI.3 of [RS3]), 

l ac (B) = wl(A,B)w+(A,B), l ac (4) = w+(A, B)w±(A, B), (4.9) 
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and the adjoints are given by 


wl(A,B) = w+(B, A) = s-lim e Its e ItA t ac (4). (4-10) 

— — t —*± 00 

Moreover, the wave operators satisfy the intertwining property (see Proposition 1 of Section XI.3 of 

[RS3]), 


wl(A,B)A= Bwl(A,B). (4.11) 

In our case, h — ho is finite rank and therefore trace-class, and we denote the wave operators w+ : = 
w+(h , ho) and the scattering matrix s := s(h,ho). 


4.2 Formulas for wave operators and the scattering matrix 

We now derive a formula for the scattering matrix. First we compute w+. We follow the methodology of 
[JKP]. Define \i = S-n and yy = Sn- 

Proposition 4.2. Let g = gi® g c ® g r e b be given. Then, 

w±g = g\ ±) ®gi ±) (4.12) 


with, 


9 ( i + \E) = gi(E) - Ji{xi, (h-E + iO) 1 g) 

O^iE) = g r (E ) - J r ( X r, (h - E + iO)" 1 ^. (4.13) 

Proof. Let any / = fi © f r e ff) ac (M = (lac (hi) © fl ac (M be given. We will compute (f,w%g). The 
computation for w* is identical. We have, 


lim (f,e ith °e- ith g) = \im(e ith e - ith °t ac (h 0 )f,g) 

t —>00 t —>00 

= (w+f,g) = < f,w\g >. 


(4.14) 


Note that, 

</, e ith °e~ ith g) = </, 5 > - i f </, e ish °ue- ish ff >d S . (4.15) 

Jo 

We require the following lemma (see Lemma 5 in Section XI.6 of [RS3]) 

Lemma 4.3. Let <p be a bounded measureable function and suppose lim t _ >00 </>(s)ds = a exists. Then 
a = lim ei0 S” e _E V(s)ds = lim £ ; 0 e _es2 </>(s)ds. 

Therefore, 


</, «’* S> = </, 5> - limi (Li(e) + L r (e)) 
£|0 


(4.16) 


where 


J l/r 


r a 

JO 


e- £8 </,e is '*» % e- is ' , 3 >d S 


J -00 

e“ e V i/r «/, e ish °5 l/r Xxi/r, e~ ish g > + </, e lsh ° X i/r)(f>i/r, e~ ish g » ds. (4.17) 
o 


Since / e hac(ho) it follows that </, e lsh °Xi/r) = 0. 
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We further compute, 


rCO 

L l/r {e)= e-e°Ji lr (f,e ish °6 l/r Xxi/ r ,e- ish g)ds 

Jo 

= Ji/r \ h/r{E) f (xi/r,e~ ls(h ~ E ~ ie) g )ds 

JR Jo 


d^l/r, ac ( E ) 


-i Ji/r \ fi/r(E)(Xl/r,(h-E-ie) 1 g)dv l/r ^ c (E) 

Jr 


(4.18) 


The interchange of the order of integration is clearly justified. Let Hi/ r (E + ie) = (xi/r, (h — E — ie)~ 1 g). 
Since lim e |o H{/ r (E + ie) =: Hi/ r (E ) exists and is finite for Lebesgue a.e. E, we have by Egoroff’s 
theorem that for any n there are measureable sets L n and R n with |IR\L n | < 1 /n, |BR\i? n | < l/n and 
Hi/ r {E + ie) —> H{/ r (E) uniformly on L n /R n . Clearly, the set 


U {fi®fr e f)ac(^o)| supp/i c L„,supp/ r c R n } (4.19) 

n>0 


is dense in t) ac (ho). Suppose that / belongs to this set. Then Hi/ r (E + ie) —* Hu r (E) uniformly on the 
support of fij r and therefore, the inequality 

\f l/r (E)\\H l/r (E) - H l/r (E + ie) | < C f \f l/r (E)\ e L\R, d^ /r , ac ), (4.20) 

which holds for a constant depending only on / and all small e implies, by dominated convergence, 

lip f \Ji/r\\Hi/r{E + ie) - H l/r (E)\dvi/ r ^ c (E) = 0. (4.21) 

£ i° Jr 

Hi/ r (E ) is bounded on the support of fu r , which implies that H l / r (E)fi/ r (E) e T 1 (1R, dvi/ r ac )- Therefore, 
the formula 


</, w*9> = ME) ( gi (E) - Mxi , (h-E- iO)- 1 ff» dv Uac {E) 

Jr 

+ f f r (E) (, g r (E) - J r ( Xr , (h-E- iO)-^)) dv r>ac ( E ), (4.22) 

Jr 

holds for the dense set of / given in (4.19). The claim follows. □ 

The following will be useful in the next subsection. 

Corollary 4.4. For a e {l,r}, 

w*h 0 S a (E) = (E5i, a + JiJ a {Xh (h- E + iO) _1 Xa> + EJ t J a (xi, (h - E + iO) -1 Xa><£a, (h 0 - E + i0) _1 <5 a >) 

© ( E5r,a + JrJaiXn (h — E + iO) 1 Xa) + EJ r J a (Xr, (h — E + iO) 1 Xa)(^a(h o — E + iO) 1 (5 a )) 

(4.23) 

and 

w*Xa(E) = (-Ji(xi, (h- E + iO) _1 Xa» © (- J r (Xr, (h- E + i0) _1 x o >). (4.24) 

Proof. (4.24) follows directly from Proposition 4.2 and (4.23) follows from Proposition 4.2 and the identity 
(Xb, (h — E + iO) -1 pA> 

= ~Mxb, (h — E + i0)-\xa> - EJ a (Xb, (h — E + iO)“ 1 Xa)<X, (h 0 - E + i0)-p o >, (4.25) 

which holds for ape {l, r}. □ 

We can now compute the scattering matrix. 
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Theorem 4.5. Let g = gi® g r e l) a c(^o) be given. Then, 


sg = gl s) ® g ( r r) 


(4.26) 


with 

gl s) (E) = 9l (E) + 2U?Ft(E)(xi, (h - E - iO^'x^E) + 2U l J r F r (E)( X i, (h-E- iO)“\ Xr ><?,(£) 

gi s \E) = g r (E ) + 2U r J l F l (E)( Xr , (h-E- iO)" 1 *,><?;(£) + 2U 2 F r (E)( Xr , (h-E-iO^xME). 

(4.27) 


There is a subtlety associated with the formulas in (4.27). For example, in the formula for g[ s '(E ), the 
function on the right g r (E ) is defined only for p r ^ c -a.e. E , and may not be defined for a set of nonzero 
pp ac measure. What we mean by the formula (4.27) is that g r = 0 outside of S ri a C . Similar statements 
hold for the formula for gi^. 

Let <v )>2 denote the standard inner product on C 2 . For any / e f) ac (ho) let f(E ) denote the vector 
(fi(E), f r (E)) e C 2 . Then for f,ge f) ac (/i 0 ) we have, 



(V(E)f(E),V(E)g(E)) 2 dE, 


where V (E) is the 2x2 diagonal matrix 


V(E) = 




(4.28) 


(4.29) 


Multiplication by the matrix V(E) is a unitary operator V : h ac (^o) ~” L 2 (R, pi(E)dE)(£>L 2 (R, p r (E)dE) 
with pij r (E) the characteristic function of 

Theorem 4.5 implies that the scattering matrix acts by multiplication by a 2 x 2 matrix on f) ac (ho). 
However, it is more convenient to consider the operator VsV -1 which is a unitary operator on Hf) ac (/io) = 
L 2 (R, pi(E)dE) ® L 2 (R, p r (E)dE) and acts as multiplication by the 2x2 matrix 


where 


s(E) 


f Su(E) Slr(E)\ 
\Sl r (E) Srr(E) J 


su(E) = 1 + 2U 2 (xi, (h-E- iO)~ 1 xi)Fi(E) 

Slr(E) = 2i JlJr{Xh (h-E- iOj-'^VWM 
Sri(E) = 2i J r Ji{xr, (h-E- iO )'- 1 X i)sjF r (E)F l (E) 
s rr (E) = 1 + 2iJ 2 < Xr , (h-E- iO )~ 1 Xr )T r (F;). 


(4.30) 


(4.31) 


This is a slight abuse of notation, and perhaps what we call s a b(E ) should really be [VsV^ 1 ^ a b(E), but 
formulas appearing later are more simple and natural with this notation. With this convention, the 
matrix s(E) is unitary for each E w.r.t. the standard inner product on C 2 . 

Since 


(Xi,(h-E- iO) x Xr) = (Xr,(h - E- iO) X xi), 


(4.32) 


the scattering matrix is symmetric for Lebesgue a.e. E e £; that is, Si r (E ) = s r i(E). From the resolvent 
identity A~ x — £? _1 = A~ X (B — A)B~ l we derive 


(Xh ( h ~ E - iO) 1 Xr) 


J?(Xh(h - E - iO) 1 xiXh, (ho ~ E- iO) 1 5 l )(xi, ( h o ~ E - iO) 1 Xr) 
1 - J 2 (br, (ho- E- iO)- 1 5 r >< Xr , (h 0 - E- iO)" 1 ^) 


(4.33) 
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The denominator is non-zero for Lebesgue a.e. E e £, as the functions (x a , (ho — E — iO) _1 Xb) are purely 
real for Lebesgue a.e. E e R. These functions are analytic outside of the spectrum of ho ['( )c (which is 
necessarily a finite set of points) and are nonvanishing. It follows from Theorems 3.17 and 5.12 of [Ja] 
that (xzj (h — E — iO) _1 xz) is non-zero for Lebesgue a.e. E e IR, and so we conclude that the scattering 
matrix is not diagonal for Lebesgue a.e. E e £. 

Proof of Theorem f.5. Let / = fi © f r ,g = gi ® g r e l) ac (ho) be given. The beginning of our proof is 
similar to the manipulations appearing in the proof of Theorem XI.42 in [RS3]. 

</, (s - t)g) = </, (w* - w* w_)g) 

= <( w+ - w-)f, w-g > 

= lim f(e ith e~ ith ° - e~' lth e lth °) f, w-g) 

£—>oo 

= lim -i f (e ish ve~ isho f,w-g)ds 
*- >0 ° J- t 

= lim—i f e- es \e ish ve- ish °f,w-g)ds (4.34) 

Jr 

The last equality is Lemma 4.3. Let us compute the integrand. We have, 

<< e ish ve~ ish °f,w_g > = J t (e- ish °f, S^e^xi, w_g) + J r (e~ ish °f, S r )(e ish Xr, w_g). (4.35) 

By the intertwining property of the wave operators, (e lsh Xi/n w -9 ) = ( w - elsh Xi/n ff) = ( w -Xi/n e~ lsh °g). 
Therefore, the formula in Proposition 4.2 gives 

</, (s - 1 )g) = limi (Hu (e) + H ri (e) + H ir (e) + H rr (e)) (4.36) 

£|0 


with, 


H ab (e ) = J a J b f e“ Es2 [ f e' sE f a (E)dv a , ac (E )1 [ f e~ isE \ X a , (h-E'~ id) Xb )g b (E')dv b ^(E') 
JlR _ JlR _ _ JlR 


ds, 

(4.37) 


Let Gab(E') = ( Xa,(h — E' — iO)xfc)- Our next manipulations are similar to the proof of Plancherel’s 
theorem in [LL] (Thm 5.3). Fubini’s theorem allows us to interchange the order of integration and 
compute, 


H r 


b (£) = J a J b f fa(E) \ f G ab (E')g b (E') [ f e-^e^-^ds 
JlR _J[R _J[R 


dv h ^ c {E r ) 


dl^a, ac(J^) 


V^JaJh J R fa(E) J G ab (E , )g b (E , ) e - 1 ^ exp du h ^(E’) 


d v a ^(E). (4.38) 


Since G ab (E')g b (E') d ^ c (E') is an i 1 (IR, d£7') function, 


f G ab (E')g b (E')e- 112 exp 
Jr 


{E 4£ J?) ) dv b ^ c (E') 2 ^G ab (E)g b (E)^p(E) (4.39) 


strongly in L 1 (R,d E) as e J, 0 (see, e.g., Thm 2.16 and Thm. 5.3 in [LL]). Let L n /R n be the set 


{E |— E^^(E) > n}. As n —* oo both \L n \ —* 0 and |i? n | — » 0. The set 

U e f)ac(/lo)|supp/i C L=,SUpp/ r C , ]| j) || ^ < GC,||/ r || 00 < GO} (4.40) 


n> 0 


is dense in i) ac (h 0 ). Fix such an /. Then the functions 

fl/r(E)^^(E) 


(4.41) 
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are bounded and so by Holder’s inequality, 

lim H ab (E) = 2J a J b f fa(E)( X a, (h-E- iO^x^F^g^dua^E). (4.42) 

£ i° Jr 

We have therefore shown that the formula 

</, (s - t)g) = </, g[ s) © gl s) ). (4.43) 

holds for the dense set of / in (4.40). The claim follows. □ 

We record here the following corollary of the unitarity of the scattering matrix; while we will not use 
it later, it is worth noting for possible future reference. 

Corollary 4.6. For Lebesgue a.e. E, 

Jr/l\(Xl/r, (h-E- iO)- l Xr/i)\ 2 F l/r (E)F r/l (E) 

= F l/r (E) Im [(J l/r (6 l/r , (h-E- iO) -1 X;/r> - l)<Xi/ r , (h - E + iO^Xl/r)] ■ (4.44) 

In particular, for vi/ r&c -a.e. E, 

Jr,l\<Xllr, (h-E- ioy'x^l 2 F r/l (E) 

= Im [(Ji/r(Si/r, (h-E- iO)“*xi/r> - l)<Xi/r, (h - E + iO) _1 Xi /r >] . (4.45) 

Proof. Unitarity of the scattering matrix implies that \su(E)\ 2 + \si r (E)\ 2 = 1. With the formulas 
established in (4.31), this implies 

J?\<Xh (h-E- iO )~ 1 Xr)\ 2 F l (E)F r (E) 

= Fi(E) (Im[<xt, (h-E- iO) _1 X;>] - JfF t (E)\{ X u (h-E- i0)-\x;>| 2 ). (4.46) 

On the other hand the identity 

(Si, (h-E- iO )~ l Xi) = —Ji(Si, (h 0 - E - iO) _ 1 <5z)<Xi, (h-E- iO) _1 Xi) (4.47) 


gives 


Im [(Ji(Si, (h — E — iO) 1 xi) ~ l)<Xb (h ~ E + iO) ^i)] 

= Im[Xi, (h-E- iO) _ 1 Xi>] - J 2 Fi(E)\(xu (h-E- iO) _ 1 X ;>| 2 


(4.48) 


which, when combined with (4.46), yields the claim. The proof of the other equality is identical. 


□ 


4.2.1 The scattering matrix for two directly coupled chains 

The formula for the scattering matrix given in [JLP] is slightly different from that derived here. We 
explain here where the difference arises. In [JLP] the XY chain consists only of a left part and right 
part and does not have a central part. To be more precise, the ho we are considering is replaced by 
ho = h r + hi where hi = h P^2(]_oo,0]) arL( I h r = h P^2(r l 00 n, and the formula given in [ ] is for s(h , ho). 

The formula obtained is that s(h, ho) acts as multiplication by the 2x2 matrix 

- (1(f) 1 { e )) < 4 ' 49 > 

where, 

ME) = 1 + 2iJg<i5i, (h-E- iO )- l h)Fi(E) 

Sir(E) = 2Uo(Jo(Si, (h-E- iO)- 1 ^) - 1 )^Fi(E)F r (E) 

sh(E) = 2Uo(Jo(S 0 , (h-E- iO)- 1 ^) - 1 )y/F r (E)Fi(E) 

s rr (E) = 1 + 2iJo<<5 0; (h-E- \0)~ l 5 0 )F r (E) (4.50) 
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with Fij r (E) = Im<£ 0 /i, ( hi/ r — E — iO) _1 d 0 /i). First note that the xi/r have been replaced by Sy .q - this 
is expected as these vectors play the same role in linking the left/right chain to the rest of the system 
through the interaction potential (i.e., these vectors appear in ‘equivalent’ places in v = h — ho and 
v = h — ho). There is also the presence of an extra term in the off-diagonal elements. This difference 
arises from the fact that when there is a central part of the chain, the \r/i are not in the cyclic subspace 
for hij r and Si/ r whereas when there is no central part, the corresponding vectors <5i ; o are precisely the 
cyclic vectors for the other part of the chain. 

The exact same methodology outlined in the proofs of Proposition 4.2 and Theorem 4.5 will yield the 
above formulas. Let us describe where the difference arises in the derivation. 

Suppose that one repeats the proof of Theorem 4.5 for h and h 0 . When one gets to (4.35) and substi¬ 
tutes for the terms (w*L(h, h 0 )S 0 /i,e~ lsh °g) (these are the terms that take the place of the (w^_Xi/rj e ~ lsh °g) 
which appear in the original derivation) one will first see a difference between the two formulas. The 
proof of Proposition 4.2 applied to h and ho yields 

w* (h, ho)So/i = [4,o/i — do((5o, (h — E + iO) 1 <5o/i)] © [4,o/i ~ Jo(Si, (h — E + iO) 1 <^o/i)] (4-51) 

where 5 a ,b is 1 if a = b and 0 otherwise. Compare this to 

w*xi/r = (h- E + i0) _1 xz/r>] © [~Jr(Xr, (h - E + iO) -1 )]. (4.52) 

The extra d 0 .o/i and <$i.o/i appearing in (4.51) lead to the slightly different formulas. 

We also note that the main results of this paper (i.e., Theorems 4.7 and 4.10) remain valid in this 
case, and the proofs carry over with little change. 


4.3 Reflectionless Jacobi matrices 


We will call a Jacobi matrix h reflectionless if the scattering matrix s(E) is off-diagonal for Lebesgue a.e. 
E e £. A priori, this definition depends on the choice of N in the formula for s(E ), but we will see in this 
section that s(E) is off-diagonal for Lebesgue a.e. E e £ for some N o s(E) is off-diagonal for Lebesgue 
a.e. E e £ for all N. 

The word reflectionless comes from the fact that the reflection coefficients |sz;(.E)| 2 and \s rr (E)\ 2 vanish 
if s(E) is off-diagonal. The reflection coefficients have the interpretation of describing the probability 
that a wave packet coming in from the left/right in the distant past is reflected and exits via the left/right 
in the distant future. There is a huge literature devoted to reflectionless Jacobi matrices and we mention 
here only [BRS], [R] and [T] and the references therein (the paper [BRS] contains a substantial list). 

There are several equivalent definitions of reflectionless appearing in the literature, and in the remain¬ 
der of this section we will discuss their relation to ours. Parts of our discussion will follow [BRS]. First 
we will require some notation. We denote the elements of the Green’s function for z e C\[R, 

gnm(z) = ( 5n , ( h - z)S m )■ (4.53) 


The limits, 


g n m{E ± iO) = lim<5 n , (h - E + ie) 1 J m ) 
£|0 


(4.54) 


exist for Lebesgue a.e. E £ IR. We denote by the operator h r« a ([ n + lj00 [) an d by h n the operator 
h [ f 2 (loc.n—il) • The Weyl m-functions are 

m n( z ) = <4+1, - z)~ l 5 n±1 ) (4.55) 

for z e C\[R, and the functions m±(E + iO) and m±(E — iO) are defined as in (4.54). In what follows we 
will sometimes adopt the shorthand g n m(E) = g nm (E + iO), g n m(E) = g nm (E — iO) to denote the limits 
whenever they exist, and the same for the mj. We have the formula [T] (see also [BRS], although their 
notation for the m functons are different than ours), 

9nU ^ J%mt(z)-m-^z)- 1 Jl_ 1 mn(z)-m+_ 1 {z)- 1 ' ^' 56 ' ) 
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If A c IR is Borel, then h is called measure theoretically reflectionless (see [13RS] ) on A if for Lebesgue 
a.e. E £ A and every n , 


Re [g n n(E + iO)] = 0. (4.57) 

Similarly, for A Borel, a Jacobi matrix h is called spectrally reflectionless on A if for Lebesgue a.e. E £ A 
and every n, 


J 2 n m+{E + iO) m“ +1 (i£ + iO) = 1. (4.58) 

Direct computation using (4.56) shows that if (4.58) holds for n and E, then (4.57) holds for n and E, 
and so, 


h is spectrally reflectionless on A ==> h is measure theoretically reflectionless on A. (4.59) 

It is well known that (see, e.g., [GKT, SY, T] and Theorem 7.4.1 of [Si] for what is probably the most 
readable proof) that 

(4.57) for E and three consecutive n => (4.58) for E and one n. (4.60) 

It is also true that [BRS] 

(4.58) for E and one n =s> (4.58) for E and all n, (4-61) 


and so h is spectrally reflectionless on A iff it is measure theoretically reflectionless on A. 
Returning to our scattering matrix, using (4.56) we find that 


s rr {E) = 1 + 2\Jhg NN {E ) ImmJ(fi) = 


Su(E) = 1 + 2iJl N _ 1 g- N - N (E) lmm_ N (E) = 
It follows that for Lebesgue a.e. E £ E riac , 


J 2 N m+(E)mN +1 ( E ) ~ 1 

J N m N+i( E ) m N( E ) - 1 

Jij V -i m Ijv(- E ) m -jv- 1 - 1 


I 


and for Lebegue a.e. E e Ej jac . 


(4.62) 


J-N-l m -N rn -N-l 1 


i ( E ) = 1 

(4.63) 

-N-l = 1 

(4.64) 


It therefore follows by the above discussion that, 


h is reflectionless oh is measure theoretically reflectionless on £ o h is spectrally reflectionless on 8. 

(4.65) 

Moreover, since, by the Kato-Rosenblum theorem, E;/ r ac does not depend on the choice of N it follows 
that whether or not h is reflectionless does not depend on the choice of N. Additionally, analoguous 
computations to those in (4.62) hold in the case of two directly coupled chains, and so the discussion 
extends to this case (in particular, our definition of reflectionless does not even depend on whether or not 
there is a central system). 


4.4 Non-equilibrium steady state 

From here on we assume that h has purely absolutely continuous spectrum. Concerning the existence of 
a non-equilibrium steady state, we have [AP], 


31 





Theorem 4.7. If h has purely absolutely continuous spectrum, then for any Ae O the limit 


<A>+ = lim lo{t\A)) 

t —>00 


(4.66) 


exists. The state w+(-) = (■)+ is called the non-equilibrium steady state (NESS) of the quantum dynamical 
system The steady state heat fluxes are 


<$;>+ = -<$ r >+ = 


i l E M E )f 


sinh(A/313/2) 


cosh(/3 T .i?/2) cosh(/3;.E/2) 


d E, 


(4.67) 


and the steady state entropy production is 


<a>+ = -/3i<$z>+ - Pr<,$r)+ = A/?<$/>+, 


(4.68) 


where A/3 = /3 r — fli. 

Proof. First suppose A e 0\ oc . The existence of the limit (4.66) is the assertion that the limit 

lim lim lim wm('Z"Lt'(^4)) (4.69) 

£-►00 M'^oo M->oo 

exists. By linearity, we may assume that A = a*(4> n )...a*((f)\)a( , if\)...a('ip n ). It follows from Proposition 
2.9 that 

lim lim (A)) = det \(e lth ipi, Te lth </,)] (4.70) 

M'—>cc M—kx) l j 

with T = (1 + e dihi+Prh r ^-i gi nce /j, 0 commutes with T and f) ac (h) = (), 

lim (e ith if,Te ith (f) = lim (e- ith °e lth ^,Te- ith °e ith (j)) 

£—>oo £—>oo 

= lim (e- ith °e ith l ac (h)if, T e - ith °e ith l a , c (h)cj ) ') 

£—>00 

= (wfip, Twfcj)) (4.71) 

since ho commutes with T. This proves the existence of (4.66) for Oi oc . An e/3 argument extends the 
result to O. 

We compute /$/)+. Our computations above have proven the formula 

w+(a*(VOa(0)) = (wlip,Tw*(t>). (4.72) 

Therefore, (2.74) implies 

<$z>+ = 2JiIm(w*h 0 Si,Tw*_xi)- (4.73) 

Let 

H ab (E) = <Xa, (h-E- iO )- 1 Xb ). (4.74) 

Direct computation using the formulas in Corollary 4.4 yields, 

2 J ; Im (w*_h 0 5 u Tw*_Xl> = J d^, ac (£) 1 ^ E (J f Im [H U (E)] - jf\H u (E)\ 2 FflE)) 

r 2 E 

- dvrME) 1 + ef3rE J 2 r Ji\H lr (. E)\ 2 F l ( E) 

-1 E ( (1 - '*"< £ " 2 > ttW - ■ < 4 - 75 » 

The formula in question then follows from the identity 1 = \su(E)\ 2 + |s; r (£')| 2 which is a consequence 
of the unitarity of the scattering matrix. The computation of /dir)-)- is similar. □ 
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If the Lebesgue measure of £, denoted \£\, is 0 then obviously there is no energy transfer between the 
left and right parts of the chain. If \£\ > 0, then since the scattering matrix is not diagonal for Lebesgue 
a.e. E G £, it follows from (4.67) that (u)+ > 0 iff fii # fi r . That is, the steady state entropy production 
is strictly positive if initially the left and right parts of the chain are at different temperature. 

In order for results similar to those of (4.7) to hold (i.e., existence of NESS, strict positivity of entropy 
production), it was necessary to take the thermodynamic limit before taking the large time limit (see 
Section 5.1 in [JOPP], and [L] for the solution to Exercise 5.1 appearing there; we repeat some of the 
discussion in [JOT ] here). More precisely, if (O', r'*, ui') is a finite dimensional quantum system, then 
the limit 

M u' 1/(4)! (4.76) 

does not exist except in trivial cases. However, the Cesaro limit 

u+(A) = ^lim w' (V‘(A)) At (4.77) 

exists for every observable A e O, and oj + is a steady state of the system. It is easy to see that 

co+(i[H',A]) =0 (4.78) 

for every observable A e O, where H' is the Hamiltonian of the system. From the above we conclude 
that it was necessary to take the thermodynamic limit of the XY chain before the large time limit. We 
also comment that the above discussion applies to infinite dimensional systems that are ‘confined’ - that 
is, H' has pure point spectrum. 

From this discussion it is reasonable to expect that we can extend Theorem 4.7 in the case that h 
has, in addition to its absolutey continuous spectrum, some pure point spectrum. Of course, the limit 
in (4.66) must be replaced with the Cesaro limit (4.77) and then the formula (4.67) holds [AJPP], The 
extension of our proof to this case is easy and requires only the Riemann-Lebesgue lemma. 

We now complete a computation which will be useful later. Let 


ko(E) 


(~PiE 0 \ 

'v 0 -B r E) ■ 


(4.79) 


Proposition 4.8. Recall the definition ofV in (4-29). Let T he a bounded operator on Uf) ac (h) that acts 
by multiplication by a 2x2 matrix 


T(E) = 


fTu(E) 

\Trl(E) 


Tlr(E)\ 

Trr(E)J ' 


(4.80) 


Then, 

tr (V~ l TVw*_ i[k, h]w-) = -J trc* (T(E)(s*(E)k 0 (E)s(E) ~ k 0 (E))) (4.81) 

Proof. Recalling that [k, h] = — fii[hi, vi\ — /3 r \h r , ty]> and using cyclicity of the trace we have , 

triV^TVwli[k, h]w-) = i ftJiiiS^hiW-V-'TVwlxi) ~ (XLW-V-'TVwlhiSi)) 

+ lfir J r ((5 r , h r W-V~ 1 TVw^_ Xr) — (XrjW-V^TVw^hrSr}) 

= iPiJi((wlh l 5 l ,V- 1 TVw*_ X i)-(w-Xi,V- 1 TVw*_h l 5 l )) 

+ i p r j r ((w*h r 6 r , V^TVwlxr) ~ {w*-Xr, V^TVwlhrSr)). (4.82) 

We would like to compute these four terms. Let H a b(E ) be as in (4.74). Note that V~ 1 TV acts on f)ac(^-o) 
as multiplication by the 2x2 matrix V~ 1 (E)T(E)V(E). Using this and the formulas in Corollary 4.4 
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to directly compute the inner products appearing above yields, 

tr (V- 1 TVwti[k,h\w-) 

= ^T n {E)E^ft i^Ji\Hu{E)\ 2 Fi{E) - 2jflm[H H (E)]) + ft2JlJ 2 \H rl {E)\ 2 F r {E)^dv^ c {E) 


+ 


iW 


Fr{E) 

Fi{E) 


E\ ft (2 J r jfH lr {E)H ll {E)F l {E) - J r JiiH lr {E )) 


+ ft (2JiJ^H rl (E)H rr (E)F r (E) + UiJ r H rl {E )) [d^ >ac (^) 


+ 


J ?;,.(£)£{&. - 2j 2 im [#„(£)]) + A2j 2 j r 2 |H ir (i;)| 2 F ; (£;)|d^, ac (£;) 

I 7 ^ 


'F,(E) 

F r {EY 


E] ft (2 J l J*H rl (E)H rr (E)F r (E) - JiJ r \H rl (E)) 


+ ft (2 J r jfH lr (E)H u (E)Fi(E) + i ftftH^E)) jdi/ riBC (£) 

= J £'|7Iz(£7) (A(|sh(^)| 2 - 1) + /3r|s;r(S)| 2 ) + 7ir(£) (, ftsu{E)sir{E ) + ftSrr{E)Srl{E)) 
+ Trr(E) (ft{\s rr {E)\ 2 - 1) + AK(£)| 2 ) + Tri{E) (fts rr (E)s rl (E) + fts lr (E)s u (E)) 


dE 
2tt ’ 


which is easily seen to be the formula in question. 

An identical computation gives, using the symmetry of the scattering matrix, 

Corollary 4.9. For T as above, 

tr(V- 1 TVw* + i[k,h]w+) = - J tr c » (T (E)(s(E)k 0 (E)s*(E) - k 0 (E))) 

4.5 The Gallavotti-Cohen functional 

The Gallavotti-Cohen functional is defined by 

GC t (a) = logw + (e - “ tEt ) 


(4.83) 

□ 

(4.84) 


(4.85) 


and describes fluctuations of the mean entropy production rate S* with respect to the NESS w + . It is the 
direct quantization of the Gallavotti-Cohen functional of classical nonequilibrium statistical mechanics 
[JPR]. Note that 

GC t (a) = lim logw s (e _Qtl;t ) = lira lim lim \oguj s = 1™ lim log uj s M (e~ at ^ tM ). 

s—>oo s— >oo M'— >oo M—> qo ’ s— > oo M— >oo ’ 

(4.86) 


Note also that the formula 

w s , M (e-“ tE ") = 


det(l + e k - s ’ M/2 e a( ' kt ’ M ~ kM ' ) e k - s ’ Ml2 ) 


(4.87) 


det(l + e kM ) 

holds. This allows us to apply the arguments in the proofs of Lemma 2.10 and Proposition 3.4 to obtain, 
Jim logw SiM (e _atE ‘ M ) = -t J dyj dutr ^1 + Sfo/ijY (4.88) 


Since k is bounded, commutes with ho and Ran(iu*) c f) ac (^-o) the strong limit 

s-lim k s = s-lim e ish e~ ish °ke ish °e~ ish = w+kwl = k+ 

s—>+GO s —> ico — 


(4.89) 
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exists. An application of dominated convergence yields, 


GC t (a) =i f dy f dutr (£ G c,t(7> u)i[k, ft]), 
Jo Jo 

with 

K GC ,t(a,u) = - (l + -u)-k- tu ) e ~k-\~ 

4.6 The entropic functionals in the large time limit 

Recall the definition of ko(E) in (4.79). Let 

K a (E) = e ko(E)/2 e a(s*(E)k 0 (E)s(E)-k 0 (E)) e k 0 (E)/2, 


(4.90) 


(4.91) 


(4.92) 


and 

K a>p (E) = ^e fco ( B )( 1 -“)/r s ( jE ) e feo(B)2a/p s ^^ e fc 0 (B)(i-a)/ P y /2 ) ( 4 .93) 


for p e]0, oo] and let 


K a00 (E) = lim K ap (E) = e (1 " a) k 0 (E)+ as (E)k 0 (E) s *(E), 

’ p —>GO ’ 

where the formula for the limit follows from Corollary 2.3 in [JOPPj. 

Since, 


, pumi= (s a (E)^E s lr (E)fi t E\ _ (s ll (E)fi l E s lr (E)fi r E\ 
1 01 JJ \s rl {E)fi r E s rr {E)fi r E) \s rl {E)^E s rr {E)fi r Ej 


(4.94) 


(4.95) 


[s(4P), fco(-E')] = 0 iff fii = fi r or E = 0 or s(E) is diagonal. We recall that s(E) is not diagonal for 
Lebesgue a.e. E e £. 

The main result concerning the large time limit of the entropic functionals is, 

Theorem 4.10. Suppose that h has purely absolutely continuous spectrum. Then the following holds: 

(i) For a e IR and p e]0, oo], 


e p+ (a) = lim -e pt (a) = 
t — >co t 


f , f det(l + K aiP (E)) \ dE 
) £ ° S Uet(l + K 0iP (E)) ) 2t r’ 

to l E S,(a) = to ]gC,(q) = e+ («) - J) log 


d E 
27T 


These functionals are identically zero iff \£\ = 0 or fii = fi r . In what follows we assume that \£\ > 0 
and fii 4= fir- 


(ii) 


The function IR 3 a i—* e Pi+ (a) is real analytic and strictly convex. Moreover, e Pj+ (0) 
—(cr) +; and 

Cp,+(o) e Pj ^.(l cr). 


-( 0 ) = 


(iii) The function IR 3 a >—» e + (a) is real-analytic and strictly convex. Moreover, it satisfies e + (0) = 0, 
e' + (0) = —(cr) + , and 


e +(0) = e" + (0) 


t^x T J 0 J t _ < cr >+)( cr _ < fT >+))+ ds | dt 


(4.96) 
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(iv) e+(l) > 0 unless h is reflectionless. If h is reflectionless then 


ia) ~ 1 


= log 


cosh((y9j(l — a) + /3 r a)E/2) cosh((/3 r (l — a) + fiia)E/2) \ d E 
cosh(PiE/2) cosh(/3 r E/2) ) 2 tt ’ 


and e+(a) = e+(l — a). 

(v) The function ]0, go] 3 p >-» e Pi+ (a) is continuous and decreasing. It is strictly decreasing for 
{0,1} unless h is reflectionless. If h is reflectionless, then e p ,+(a) does not depend on p and is 
equal to e+{a). 

Remark. The time-reversal invariance of our system together with Corollary 4.9 implies that we can switch 
s(E) and s*(E) whenever they appear in the equations (4.92)-(4.94) when inserted into the formulas in 
part (i). This is can also be seen directly by conjugation because s(E) = s*(E). 

Proof. We begin with part (i). Note the existence of the strong limits 


IC+(a) = s-lim/CES t{ot,u) = s-lim/Ccc t(oi,u) = — (1 + e 

t —>00 ’ t —XX) ’ V 

and, for p e]0, oof, 


-a(k + -k-) e ~k 


C ") 


-l 


(4.97) 


/C Pi+ (a) = s-lim/C Pi t(a, u) 


= _ P -( 1 -“) fe +/p 
2 


^1 + ^ 1 ~ a ^ k +/P e 2o ‘ k -/P e ^ 1 ~ a '> k +/p^- p / 2 '\ 1 e ( 1_a ) fe +/p + h.c. 


and 


JCco + (a) = s-lim/Coo,t(Q!,u) = (t + e (1 a)k+ ak ~\ 
t-> 00 V / 


(4.98) 


(4.99) 


This follows immediately from the definitions of the and k ±, where ff stands for ES, GC or p e]0, oo]. 
By (3.21) and (4.90), 

i nCt s* 1 s>Ot n\ 

-e#, t (a) = dy tr (/C #i t(y, u)i[k, h]) = V dy du(6 k , £#,t(y, u)\[k, h]S k ), (4.100) 
1 Jo Jo keA^° 

where ecc,t(a) := GC t (a) and A = [—IV— 2, —TV +1] u [TV— 1, TV+ 2] c Z. The aforementioned existence 
of the strong limits and dominated convergence implies, 

lim ie #jf (a) = ^ f <4,/C #i+ (y)i[fc, h]4> 

°° 1 fce_4 Jo 

= f dytr(/C #i+ (y)i[fc,h]). (4.101) 

Jo 

In particular, 

e + (a) = lim -ES t (a) = lim -GC t (a) = f tr(/C + (y)i[fc, Ti])dy. (4.102) 

t —>00 t t-KX) t Jg 

By the dehnition of the scattering matrix and Proposition 4.8, 
tr(/C+(y)i[fc, h\) = - tr((l + ^ 

= J tr ((1 + e-^^* {E)ko{E)siE) - ko{E)) e- koiE) )- 1 {s*{E)k 0 {E)s{E) - fc 0 (E))) (4.103) 
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At this point, we can essentially reverse the computation completed in Lemma 2.10. By Corollary 2.12, 


tr f(l + e -^*(E)ko(E)s(E)-k 0 (E)) e -k 0 (E)^ 1 ( s * (£) (E)s(E) - k 0 (E))\ 

= — tr log (1 + e ko(E)/2 e j( s *(E)k 0 (E)s(E)-k 0 (E)) e k 0 (E)/2^ _ 

For matrices A and B the equality 1 + e A e B = e -B//2 (l + e B / 2 e A e B / 2 )e B / 2 implies 


(1 + e -y(s*(E)ko(E) s (E)-ko(E)) e -ko(E)yl 

< exp 

max \(3 a E \ 



Ee £, ae{Z,r} 


and so, Fubini’s theorem yields (recall that £ is bounded), 


(4.104) 

(4.105) 


(4.106) 


e+(a) = f tr(£+( 7 )i[fc,/i])d 7 
Jo 

= J — tr log (l + e ko(E)/2 e^( s *( E)ko{E)s{E) - ko{E)) e ko{E)/2 ^ dy 

det(l + e ko(E)/2 ea (s*(E)k 0 (E)s(E)-k 0 (E)) e k 0 (E)/2^ dE 


d E 
27r 


1 


= log- 


det(l + e k °( E )) 
which is the formula in question. For p e]0, oo] 


27r 


(4.107) 


tr(£p 1 + ( 7 )i[fc,/i]) = “ J tr(7;, P (S)(«*(£0*o(^a(£0 - k 0 (E))) ^ (4.108) 

where, for p < oo, 

7 ^ p (E) = ie _< ' 1 _ 7 l s * fcoS ^' E ^ p + (e^ 1 _ 7 l s * feoS l■ E l/ p e 27 fco ^ E l/ ^ ’el 1 _ 7 l s * fcoS l £: l/ ^ ’) _p/,2 ^ 1 e ( 1 ~'r) s * k os(E)/p 

+ h.c., (4.109) 


and 


7^,00(-E) = ^l + e" (1 " 7)s * (i5)fco(B)s( ' B) “ 7fco(B) ) \ (4.110) 

and s*kos(E ) = s*(E)ko(E)s(E). By Corollary 2.12 for p e]0, oo[, 

ti (Ty, p {E)(k 0 (E) - s*{E)k 0 (E)s{E))) = 4-trlog (l + ^-As*k oS (E)/ Pe 2 jk 0 (E) e (i-^ s *kos(E)/ P y /2 ^ 

= A- tr log (l+ ^ e ( 1 - 7 )fco(E)/p s(F;)e 27 feo(E) s*(F;)e (1 - 7)feo(E)/p ) P/2 ^ 
= ^-trlog(l + K 7:P (E)) , (4.111) 

where the second equality is just the fact that s(E) is unitary. We have also, 


-ti(Ty,o,(E)(s*(E)k 0 (E)s(E) - k 0 {E))) = J-logtr (l + e -^-'y)s‘(E)k 0 (E)s(E)- 1 k 0 (E)^ 

= A. log tr (l + e A-'y)ko(E)+'ys(E)k 0 (E)s*(E)^ 

= Ai 0 gtr(l + AT 7 , 00 (E)). (4.112) 
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It is easy to see that, for p < oo, 


||7^, p || < exp 


|1 — 7b 1 max \PaE\ 

EeS, ae{Z,r} 


(4.113) 


and 


II ,00 II ^ 1; 


(4.114) 


and so Fubini’s theorem yields for every p, 


e p,+ 



tr (£ P ,+ (7)i[Mhl)d7 



d 

dy 


tr log K 1 ' P (E)&'y 


d E 
2i r 


f det(l + K a}P (E)) d E 

J £ ° g det(l + Kq(E)) ~2tt 

(4.115) 


Note that if G(a ) is a differentiable matrix-valued function taking values in invertible matrices, then, 

4-G~ V) = —G~ 1 (a)G'(a)G~ 1 (a). (4.116) 

da 

This and dominated convergence yield, after some simple algebra, 

e" (a) = J -d- tr (e ko{E) [e k °^ E) + e -<*(s*(E)k 0 (E) s (E)-k 0 (E))^ _1 ( 8 *(E)ko(E)s(E) - k 0 (E)) \ ^ 

= J tr ((e“ A(B) + e- k ^ E) )- 1 A(E)(e~ aAiE) + A{E)} ^ (4.117) 

where A(E) = s*(E)ko(E)s(E) — ko{E). The integrand is non-negative for every E and vanishes 
iff A(E) = 0 o [s(E),ko(E)] = 0. Since e+(0) = 0, it follows that e+(a ) vanishes identically iff 
[s(.E), ko(E)] = 0 for Lebesgue a.e. E e £ or \£\ = 0. 

An explicit computation shows that [s(.E), ko(E)] = 0 iff Pi = f3 r or s(E) is diagonal. Since s(E) is 
not diagonal for Lebesgue a.e. E e £ it follows that e + (a) vanishes identically iff |£| = 0 or /3; = /3 r . 
Similarly, with B{E) = s*(E)ko(E)s(E), 

e 2,+ ( a ) = - J ^ tr ( [(1 + e-^- a)B(E) e- a k 0 {E))~ l + (1 + e -^ 0 (E) e (i- a )B(E)yij 

) d F 
4^ 

= J tr (l e A~ a ) B ( E ) + e - ako( ' E) )- 1 (B(E) - k 0 (E))(e ako{E) + e -G-«)B(E)yi 
x (B(E) - ko(E)) + h.c.) (4.118) 


The integrand is non-negative and vanishes iff [s(£'), ko(E)] = 0 and so as in the case of e+(a), we 
conclude that e 2 ,+ (a) vanishes iff |£| = 0 or /?/ = /3 r . 

If e 2 ,+(a) vanishes identically, then [s(.E), fco(-E)] = 0 and it follows that e Pj+ (a) vanishes identically. 
On the other hand, if e 2 ,+ (a;) doesn’t vanish identically then the above computation shows it is a strictly 
convex function. In (ii) we will see that e 2 ,+ (l) = e2,+(0) = 0 and so e 2 ,+ (a) < 0 for a e]0,l[ and 
e 2 ,+ (a) > 0 for a $ [0,1], In (v) we will see that the function ]0, oo] 3 p —* e Pi+ (a) is decreasing, and so 
it follows that e Pj+ (a) does not vanish identically if e 2 ,+ (a) does not vanish identically. This completes 
the proof of part (i). 

We now prove parts (ii) and (iii). The analyticity of the entropic pressure functionals follows directly 
from the analyticity of the functions 


(5 k ,)C# !+ (a)i[k,h]5 k ). 


(4.119) 
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The same arguments which led to the analyticity of the finite time entropic functionals in Proposition 3.4 
yields that the above functions are analytic (note that here the argument is easier as there is no integral 
and no t, so the Fj’s will be taken to be functions only of a). 

Convexity and the symmetry e Pi+ (a) = e Pi+ (l — a) follow from the fact that these properties are 
satisfied by the finite time functionals, of which these are pointwise limits. The fact that e Pi+ (0) = 0 
follows from the formulas in (i), or that the finite time functionals satisfy the same equality. The formula 
(4.117) shows that e + (a) is strictly convex unless either \£\ = 0 or [s(£), fco(S)] = 0 for Lebesgue a.e. 
E e £ iff ft = (3 r . 

As limits of convex functions, the e Pj +(a) are convex. Their second derivatives are analytic and 
therefore are either identically 0 or have a set of isolated 0’s. If they vanish identically then the e Pj +(a) 
are linear; however, since e Pi+ (l) = e Pj+ (0) = 0, they then must vanish identically if they are linear. 
Hence, the e Pi +(a) are strictly convex if they do not vanish identically. 

Since /C Pi +(0) = (1 + e~ k+ )~ 1 for p e]0, go] and /C+(0) = —(1 + e~ k ~)~ 1 , we have by the first line of 
(4.107) 

e' + (0) = — tr ((1 + e _fe_ )^ 1 i[fc, h]) = — tr ((1 + e~ k )~ 1 w*L i[fc, h]w-) = —(cr) + , (4.120) 

and by (4.115), 

e' p + = tr ((1 + e _fc_ ) _1 i[fc, h]) = tr ((1 + e _fe ) _1 w*i[fc, h]w+) = — tr ((1 + e~ k )~ 1 w^_ i[fc, h]w-) = —(&)+, 

(4.121) 


where the second to last equality follows from the time reversal. 

The application of Lemma A.l in the proof of Proposition 3.4 yields that there is an e > 0 so that the 
functions ES t(a)/t and e 2 ,t{oi)/t have analytic extensions to the disc D{ 0,e) and are uniformly bounded 
in t > 0 on this disc. The same arguments apply to the functions GC t(a)/t. The Vitali convergence 
theorem implies, 


s 2 +( 0 ) = lim -e 
• +v ' t-00 t 


( 0 ), 


el(0) = lim yES"(0) = lim JgC"(0). 
t— >00 t t— >00 t 


(4.122) 


By Proposition 3.4, e^'^O) = ES"(0) for each t and so e2 + (0) = e" (0). By the Vitali convergence 
theorem, 


i2 

GC"(0) = lim lim f- loguj sM (e~ at ^) 
s —>00 M—>oo da 2 \ / 


(A = 0 


Direct computation yields 
d 2 


-j -2 log W Sj M (e QtE ‘ M ) = w Sj M - (w Sj M (tE^)) 2 

(&u,m')^s,m (&v,M')dudv, 

Jo Jo 

and so by dominated convergence, 

7 GC"( 0 ) = \ f f lim lim (ws^K.m^.m) - w s ,mK,m)w s ,m(<t»,m)) d«d« 
t t Jo Jo s—>G0 M—>oo 

= 7 I I (a u a v )+ - (a u ) + (a v ) + dudv 

t Jo Jo 

n (<7u- v cr)+ - <cr)+dwdw 
) 

J t «^cr>+ - <( 7> 2 ) ^1 - \~Pj ds 

J <(cr„-<0-> + »(cr-<f7> + )> + ds. 


(4.123) 


(4.124) 


(4.125) 
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In the third and fifth lines we have used w + o r 4 = u> + , and the fourth line follows from a change of 
variable. Integration by parts yields, 


1gc"(0) 



-<c r )+)(o--<o')+))+dM 


ds, 


(4.126) 


which completes the proofs of (ii) and (iii). 

We turn to the proof of (iv). We will compare e+(l) and eoo(l). For any 2x2 matrix, det(l + A) = 
1 + tr(.A) + det (A). We have also, 


det ( K atCO (E )) = exp [tr ((1 - a)k 0 (E) + as(E)k 0 (E)s *(£'))] 

= exp[tr (fc 0 (.E))] 

= det (. K a (E )), 

The equality 

tr ^ e k o(E)+ a (s(E)k 0 (E)s*(E)-k 0 (E))^ = tr ^ e k 0 (E)+ a (s*(E)k 0 (E) s (E)k 0 (E)-k 0 (E))^ 

follows by conjugation and so, 

det (1 + K a (E)) - det (1 + K a>o0 (E)) = tr f e ko(E) e aA(E)\ _ tr ^ e k 0 (E)+aA(E)\ ^ 


(4.127) 

(4.128) 


(4.129) 


with A(E) as before. By the Golden-Thompson inequality (see Corollary 2.3 and Exercise 2.8 in [JOPP]; 
see also [L]), the RHS is strictly greater than 0 unless s*(E)ko(E)s(E) and ko(E) commute. It is easy 
to see that, if A /3 r , this can happen only if s(E ) is diagonal or off-diagonal or E = 0. For example, 
one can compute, 

[k 0 (E)s*(E)ko(E)s(E)] h2 = E 2 (tfs u (E)s lr (E) + ^ r s rr {E)s rl {E)) (4.130) 

and 


[s*{E)k 0 (E)s(E)k 0 (E)] ly2 = E 2 W r su{E)sir{E) + ffs rr (E)s rl (E)). (4.131) 


By unitarity of the scattering matrix, su(E)si r (E ) + s rr (E)s r i{E) = 0. It follows that if ko{E ) and 
s*(E)k 0 (E)s(E ) commute, then either E = 0 or /3/ = /3 r or su(E)si r (E) = 0. Since the scattering 
matrix is unitary, this can happen only if s(E) is diagonal or off-diagonal. Since s(E) is not diagonal for 
Lebesgue a.e. E e £, it follows that s*(E)ko(E)s(E) and kg(E) commute for Lebesgue a.e. E e £ iff 
s(E) is off-diagonal for Lebesgue a.e. E e £, i.e., iff h is reflectionless. In summary, we have proven that 


e+(l) > eoo, + (l) = 0 (4.132) 

and equality holds iff h is reflectionless. 

If h is reflectionless, then, 

s*(E)k 0 (E)s(E) = E , (4.133) 

and so, 

det (1 + K a (E)) = (1 + e -/3rBa-( 1 -«)/3 i B)( 1 + e -fhE a -(l- a )0 r E ) 

_ 4cosh(G9,(l -a) + /3 r a)E/2) cosh((/3 r (l - a) + fta)E/2) fA io n 

e (f3i+p r )E /2 ’ (. 4 . 1 ^ 4 ; 

from which the formula in (iv) follows. 

Only (v) remains. We have, 

det (1 + K a , p {E)) - det(l + K a<q {E)) = tr {K a , p {E)) - tr {K a<q (E)) . (4.135) 

The Araki-Lieb-Thirring inequality (see Theorem 2.2 and Exercise 2.8 in [JOPP]; see also [L]) implies that 
the RHS is strictly positive for p < q unless s*(E)ko(E)s(E) and ko(E) commute iff h is reflectionless. 
If h is not reflectionless, it follows that ]0,oo] 3 p —* e Pj +(a) is strictly decreasing. If h is reflectionless 
then it is trivial to check that all these functionals are identical and equal e+(a). □ 
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4.7 Large deviations 

We discuss here some consequences of Theorem 4.10. In this section we assume that \£\ > 0 and /3; ^ /3 r . 
Recall that, 

e 2lt (a) = FCSda) = log f e^dF^), (4.136) 

Jr 

where P t is the FCS measure of the extended XY chain of Proposition 3.4. In the GNS representations 
for O and the states w and w+ (see, e.g., [BR1]) denote the spectral measure for X* and w and w + by 
P E s,t and P GG ,t, respectively. Then, 

ES t (a) = log f e- a ^dP ES ,*(</»), GC t (a) = log f e -“^dP GC ,*(</>)• (4.137) 

Jr Jr 

The large deviation rate functions are given by, 

Jfcs+($) = - inf (a6 + e 2 ,+ (a)) 

aeR 

I+(6) = — inf (a6 + e+(a)). (4.138) 


We have, 

Lemma 4.11. (i) The large deviation rate functions are non-negative, real-analytic, strictly convex 

and vanish at the single point 9 = (</>)+. 

(ii) The two rate functions are different unless h is reflectionless. 

Proof. The results follow from basic properties of Fenchel-Legendre transforms and the facts that e' + (0) = 
e 2 +(0) = — (4>y+ and e + (l) > 0 = e 2j +(l) if h is not reflectionless. □ 

The symmetry e 2 ,+ (a) = e 2 ,+(l — a) implies , 

Ifcs+(6) = Jfcs+(“$) + (4.139) 


If I+(6) satisfies this relation, then the symmetry e+(a) = e+(l — a) must hold and so h is reflectionless. 
The consequences of Theorem 4.10 are, 


Corollary 4.12. Suppose that h has purely absolutely continuous spectrum. 

(i) The Large Deviation Principle holds: for any open set OcR, 

lim -j- log P E s 1 t(0) = lim \ log P G c,t(0) = - inf I+(9), 

t->oo t t —>oo t OeO 

lim j log P FG s,t(0) = - inf / FC s+(6 ) )- 
t-KX) t OeO 


(ii) The Central Limit Theorem holds: for any Borel set B c R ; let B t = {(f> \ — e B}. Then 

1 


lim P E s ,t(B t ) = lim P GG ,t(S t ) = lim P F CS,t(Rt) = ,-— 

t ->00 *-► 00 ’ t — 00 ’ W2nD. 


: f e-^^dcf, 

Jb 


where the variance is D + = e" (0). 

Proof. Since e + (a ) and e 2 ,+(o:) are, in particular, differentiable, part (i) follows from the Gartner-Ellis 
theorem (see, e.g., Appendix A.2 in [JOPP]). In the proof of Theorem 4.10 we showed that there is an 
e > 0 so that the functions e 2 ,t(a)/t, ESt{a)/t and GCt(a)/t have uniformly bounded extensions to the 
disc D( 0,e). This, with the fact that e' + (0) = e' 2 + (0) = — (<j)+, implies that part (ii) follows from Bryc’s 
lemma (see Appendix A.4 in [JOPP]). □ 
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5 Conclusions 


The differences between quantum and classical mechanics is due to the non-commutative structure of 
quantum mechanics. In quantum statistical mechanics this results in the emergence of novel entropic 
functionals. In our notation, the functionals e Pt t(a) and ESt(a) all result via different quantizations of 
the classical entropic functional. Note that the entropic pressure functionals all have the a <-» 1 — a 
symmetry which has played a key role in the recent developments in classical non-equilibrium statistical 
mechanics, while the direct quantizations of the Evans-Searles and Gallavotti-Cohen functionals do not. 

Theorem 4.10 shows that the Gallavotti-Cohen and Evans-Searles functionals are identical. This is 
not surprising, as the same phenomenon will holds in classical statistical mechanics - this is known as the 
principle of regular entropic fluctuations - see [JPR]. The fact that the entropic pressure functionals all 
equal the Evans-Searles functional if the Jacobi matrix is reflectionless is remarkable. Recall that ei jt (a) 
was linked with the full counting statistics of a repeated quantum measurement protocol in Section 
2.3.3. This measurement protocol is purely of quantum origin, and has no classical analogue. That this 
functional should equal the one arising from classical mechanics is surprising. 

Recall that our hypotheses required that the Jacobi matrix h have purely absolutely continuous 
spectrum. The case that is most important for quantum mechanics is the Schrodinger case, in which all 
the J n ' s are all equal (recall that the J„’s are the nearnest neighbour coupling of the spins). What is 
remarkable is that all known examples of this type of Jacobi matrix with purely absolutely continuous 
spectrum are reflectionless. For example, if J n = 1 and v n = 0 for every n, then the resulting operator 
is called the Laplacian, denoted A. The discrete Fourier transform identifies A with the operator of 
multiplication by 2 cos 9 on L 2 ([—7r, 7r], d0), and so A has purely absolutely continuous spectrum. Explicit 
computation shows that A is reflectonless (see, e.g., [L]). 

To summarize, it is currently an open problem if there exists a Jacobi matrix with J n = J ¥= 0 with 
purely absolutely continuous spectrum that is not reflectionless, and it would seem that this problem 
is relevent to the theory of entropic fluctuations. It is believed that such an example exists [JL ]. In 
the non-Schrodinger case it is easy to construct operators that have both purely absolutely continuous 
spectrum that are not reflectionless - see [VY]. 

In the remarks following the proof of Theorem 4.7, we noted that the existence of a NESS could be 
extended to the case where h has some pure point spectrum if the steady state is replaced by a Cesaro 
sum. It is currently open as to whether a result of this sort can be extended to Theorem 4.10. That is, 
do the large time limits of the entropic functionals exist if h has some point spectrum? 

There are few models in quantum mechanics for which the existence of the large time limit entropic 
functionals can be proven. The XY chain is a useful model because the functionals can be computed 
in closed form in terms of the scattering data and their properties can be examined. It remains to be 
seen that if the identification of XY chains with identical entropic functionals with reflectionless Jacobi 
matrices can be extended beyond exactly solvable models like the XY chain in the future. 

A Analyticity lemma 

We record here and prove a lemma that was used to obtain the analyticity of the entropic functionals. 

Lemma A. 1. Let f{z) be a function analytic on the half-plane Re z > 0. Let Ffla, s, t) : IR x [0,1] x [R —» 
O, j = 1,2,3 take values in bounded operators on a Hilbert space. Suppose that there is a strictly positive 
continuous function a(a) s.t. 


Fi(a, s, t) > a(a)t (A.l) 

for every s and t. Suppose that for fixed s,t the function a —» Fj(a,s,t) is analytic in the following 
uniform sense: for each ao e IR, there is an e > 0 so that the series 

00 

Fj(a,s,t) = A n j(s,t)(a- a 0 ) n (A.2) 

n=0 
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is absolutely convergent for |a — ao| < e. Suppose furthermore that each function s —» A n j(s,t) is 
continuous, and that there are numbers a n j s.t. ||u4 ni j(s, f)|| < a n j for every s,t and j and the function 

CO 

H A°) = X a n,j\a - cto\ n (A.3) 

n=0 


is finite for \a — «o| < £• 

Then for any -0 and f>, the functions 

9t(oi) = f ('ip,F 2 (a,s,t)f(F 1 (a,s,t))F r (a,s,t)(/>)ds (A.4) 

Jo 

are real-analytic in a. Furthermore, for each ao e IR, there is an e > 0 so that each gt(a) is analytic in 
the disc D(ao,e) and 


sup \gt(a)\ < co. (A.5) 

aeD(ao,e) 
t 

Proof. Fix ao e R and let e > 0 be as in the hypotheses on the Fj and H. WLOG we can suppose that 
a n ,i = <in ,2 = On ,3 for every n and we drop the subscripts and write a n and H for a n j and Hj. Let, 

M = sup H(ct) < go, a = inf a(a) > 0. (A.6) 

aeD(a o,e/ 2 ) ae[a 0 -E,a 0 +£] 


Choose X > M + a. Define, 




H(a) = X — a + {H(a) — do) > 0. 

(A.7) 

It follows that, 




||Fi(a, s, t) - A | < H(a), 

(A.8) 


for a e D(ao, e). Since H(ao) = X, there is an 0 < e' < e/2 and a S > 0 s.t. 

H{a) < X — 6, (A.9) 

for a e D(ao,e'). Since f(z) is analytic in the half-plane Rez > 0, it follows that the series 


CO 

f(z) = 2 b n {z - A)" (A.10) 

n=0 

converges absolutely for \z — X\ < X. Note that s’ does not depend on s or t. It follows that f{F\{a, s, t)) 
is analytic for a e D(ao,e'). 

The formal power series obtained by expanding each of , s, t), f(Fi(a, s, t), Fs(a, s , t), sandwiching 

them between then ip and the <f> and integrating over s in fact converges absolutely in the disc D(ao,s') 
because, 


f 1 m \\n 

00 

X l|A„, 2 (s,t)|| | a-a 0 \ n 

00 

X IK, 3 (M)II 

a - a 0 \ n 

Jo 

_n =0 

n=0 


00 

00 


m 

x X! I 6 ™! 

||A 0 ,i(s, t) - A|[ + X l|An,i(s,t)|| \a - a 0 | n 

d s 

o 

II 

n= 1 

00 

00 



< Ill’ll H\\ H(a)H(a) J] \b m \H(a) m < \\iP\\ \\<P\\ M 2 J] \b m \(X - S) m < co. (A.ll) 

m =0 m =0 

This yields the claim. □ 
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